WEIGHTED BERGMAN SPACES: SHIFT-INVARIANT 
SUBSPACES AND INPUT/STATE/OUTPUT LINEAR SYSTEMS 
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Abstract. It is well known that subspaces of the Hardy space over the unit 
disk which are invariant under the backward shift occur as the image of an 
observability operator associated with a discrete-time linear system with sta- 
ble state-dynamics, as well as the functional-model space for a Hilbert space 
contraction operator, while forward shift-invariant subspaces have a represen- 
tation in terms of an inner function. We discuss several variants of these 
statements in the context of weighted Bergman spaces on the unit disk. 



1. Introduction 
Let H 2 be the standard Hardy space over the open unit disk 

oo oo 

ff 2 = {/(*) = £/»*"= ]T|/„| 2 <oo} 

and we let S: f(z) >— > zf(z) be the shift operator on H 2 . The classical theorem of 
Beurling [TT] asserts that any ^-invariant subspace A4 of H 2 has the form M = 
9 ■ H 2 where 9 is an inner function (analytic on the disk with boundary-value 
function on the unit circle having modulus 1 almost everywhere). The result was 
extended to the vector-valued case by Lax [32] (for the finite-dimensional case) and 
by Halmos [IB] (for the general case). The main thrust of Beurling's paper was the 
development of a theory of inner-outer factorization for H 2 and H °°-functions; this 
theory was then used to arrive at the famous characterization of invariant subspaces 
M = 9-H 2 . 

There are now a number of more operator-theoretic and/or system-theoretic 
proofs which better handle the vector-valued case. For y a Hilbert space, we write 
H 2 (y) for the Hardy space H 2 (y) = H 2 ®y of y~ valued functions and we let M be 
a shift-invariant subspace of H 2 (y). We single out four distinct approaches toward 
what we shall call simply the Beurling-Lax theorem for the vector-valued case. 

(1) Find a Hilbert space X (playing the role of a state space from the system 
theory point of view) and construct operators C : X — > y and A : X — > X so that 
[q] is isometric (A* A + C*C = Ix) and so that M 1 ^ has the representation 

M 1 - = {C{I - zA)- 1 x: x G X}. (1.1) 

One such choice is X = M.^ with C = E\ M ± (where E: f M> /(0) is evaluation at 
zero and A — S*\ M ±). We note that H 2 (y) is a reproducing kernel Hilbert space 
with reproducing kernel given by the Szego kernel Then M 1 - C H 2 (y), as a 
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subspace of H 2 (y), is also a reproducing kernel Hilbcrt space in its own right. One 
can identify the reproducing kernel for .M^ explicitly as 

M ± =U{K c ,a) where K c ,a = C(I - zA)~ l {I - (A*)' 1 ^. 

It then follows that M. has Beurling-Lax-Halmos representation M = • H 2 (U) 
where is any solution of the kernel factorization problem 

0(^)0(0* h K , h 9 , 

— ; 7— = - j-Kc,a{z,Q (1.2) 

1 — zQ 1 — zQ 

from which it follows that the reproducing kernel kj^{z, () for the subspace M. is 

0(z)0(Q* n 
km(z,C) = — =— • (1-3) 

1 - < 

One way to construct such a is as 

e(z) = D + zC{I - zAY x B (1.4) 
where [ £ ] solves the Cholesky factorization problem 



[B* D*} 



"I 


0" 




A 





/ 




C 



[A* C*] 



We refer to this approach as the de Branges-Rovnyak-Potapov approach, as the 
first part (the identification of M. as W(K c ,a ) is prominent in the work of de 
Branges-Rovnyak H3] while the second step (factorization of the kernel (|1.2[) ) 
is prominent in the Potapov-school approach to interpolation (see |17|). 

(2) Set A = PmS*\m and C = Ps\m> M = M Q SM and consider the operator 
O c ,a :McH 2 (y)^ H 2 {U) given by 

O c ,a - C(I - zA)- x x. 

Then the inner representer for M can be obtained as 

M e = (O c ,a)* : H 2 {U) -> H 2 (y). 

This approach works more generally for S*-invariant subspaces of H 2 (y) contrac- 
tively included in H 2 (y). We refer to this approach as the Rosenblum-Rovnyak 
approach, as this is amounts to the proof of the Beurling-Lax Theorem in the book 

ESI. 

(3) Note that M decomposes as M = fc>o (S k M 6 S k+1 M). It turns out 
that the reproducing kernel for the subspace S k A4 S k+1 Ai has the form 

ks>MeswM(z, = z k ( k e(z)e(cy (1.5) 

with Q(z) independent of k. Thus the operator of multiplication by z k Q(z) maps 
the coefficient space U isometrically onto S k M. Q S k+1 M. and we recover the kernel 
function for M to be 

00 

W*,C) = 5> fe f 0(2)0(0* 

in agreement with (ll.3[) . It now follows that is a Beurling-Lax representer for 
M.. We refer to this as the iterated Halmos wandering- subspace approach since 
the construction of the inner representer via looking at the wandering subspace 
M. SM is the main idea of the proof for the abstract shift-operator setting in 
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(4) Construct 9 so that M&U = M SM, or k M eSM = Q(z)Q(()* (i.e., 
is the same as in Approach 3). In this shift-operator context, then the mul- 
tiplication operator Mq extends to be an isometry from H 2 (U) into H 2 (y) with 

MqH 2 (U) = M and once again k M (z,C) = ®(z) (7^) ©(C)*; so once again 
is the Beurling-Lax representer for M.. We refer to this as the uniterated Halmos 
wandering subspace approach. 

In all four approaches, one ultimately arrives at the realization formula (jl.4p for 
the inner function 9 where the system matrix (also called operator colligation) U = 
[ £■ p ] is unitary. This has the interpretation that the associated input/state/output 
linear system 

J x(n + 1) = Ax(n) + Bu(n) , . 

U ' \ y(n) = Cx(n)+Du(n) { J 

is conservative: the energy stored by the state at time n (\\x(n + 1)|| 2 — ||cc(n)|| 2 ) 
is exactly compensated by the net energy put into the system from the outside 
environment (||u(n)|| 2 — ||y(n)|| 2 ), with a similar property for the adjoint system. 
We note that application of the Z-transform 

00 

{x(n)} ne % + i-> x(z) := ^ x{n)z n 

converts the system equations (jl.6|) to 

z^ 1 (x(z) - x(0)) = Ax(z) + Bu(z) , . 

y{z) = Cx(z) + Du(z) 1 '> 

which can then be solved for x{z) and y(z): 

x{z) = (I — zA)~ l x{0) + z(I — zA)~ l Bu{z) 
y{z) = (Oo,ax(0)) (z) + e v (z) ■ u(z) 

where 

O c ,a- x^C{I -zAY x x (1.8) 
is the observability operator for the system Eu (generating the Z-transform y(z) of 
the output signal from the initial state x(0) = x when the input signal is taken to 
be zero), and where 

9u(z) = D + zC{I - zA)- 1 B (1.9) 
is the transfer function of the system Su (having the property that multiplication 
by 0u on the Z-transform u{z) of the input signal generates the Z-transform y{z) 
of the output when the initial state x(0) is taken to be zero). We note that the 
observability operator Oq,a appears in the representation for A4 and that 
the formula (11.41) exhibits the Beurling-Lax representer for Ai as the transfer func- 
tion for the conservative i/s/o linear system Su- Many of these ideas connecting 
operator-valued H 00 -functions with systems theory ideas can be found already in 
the survey paper of Ball-Cohen [10] . 

It was only much later that researchers began looking for analogues of Beurling- 
Lax representations for shift-invariant subspaces of weighted Bergman spaces (see 
[T^l [T5] and the references there for the scalar- valued case) . For a Hilbert space y 
and an integer n > 1, we denote by A n (y) the reproducing kernel Hilbert space 
with reproducing kernel 

K(z,0 = j: 1 7 ly (1-10) 

(1 - z() n 
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which is clearly positive on B x ID). The space A n (y) can be alternatively charac- 
terized as the Hilbcrt space of J^-valued functions / analytic in the open unit disk 
D and with finite norm ||/||.A„(;y) : 



My) = { /(z) = E-^' : H/ll2uoo : =E^ • ll/illy < 00 \ ( L11 ) 



We let S n denote the Bergman shift operator of weight-index n acting on A n (y) 
given by multiplication by the coordinate functions: S n : f(z) h-> zf(z) for / G 
A n {y). The space A\{y) is therefore the standard Hardy space H 2 (y) of the unit 
disk with S\ — S equal to the standard Hardy-space shift operator, and for n > 2 
the space A n (y) is the standard weighted Bergman space. 

Even for the case y — C, one can no longer use Blaschke products to collect the 
zeros of a Bergman-space function since the zero set of a Bergman space function 
need not satisfy the Blaschke growth condition. A major advance came with the 
work of Hedenmalm |20) who constructed Bergman- inner functions, i.e., a function 
9 so that 9C = A2 S2A2 as the solution of an extremal problem. Such a 6 
has the contractive divisor property \\9 ~ 1 f\\A 2 ^ 11/11-42) thereby improving earlier 
results of Horowitz [3T| . Since the work of Apostol-Bercovici-Foias-Pearcy [5] it is 
known that shift-invariant subspaces M. C A2 can have arbitrary index indA^ := 
dim (Ai © S2M). Nevertheless, the seminal work of Aleman-Richter-Sundberg [2 
with later extensions by Shimorin [31] showed that in all cases we recover M. as 
Ai = spaE k>0 S^(Ai S n A4) (a partial analogue of the representation of M. in 
approach 4 above), at least for the cases n = 2,3. More recently, Olofsson [2"T1 
1281 129) initiated the study of operator- valued Bergman- inner functions for standard 
weighted Bergman spaces as an object of independent interest. In particular he 
suggested a certain time-invariant input/state/output linear system of higher order 
(past values of the states and of the inputs enter the state-update equations rather 
than just the current state and input values) as the time-domain explanation for the 
input-output map corresponding to multiplication by the Bergman inner function. 

The purpose of the present paper is to further enhance the theory of Beurling-Lax 
representations and to sharpen the connections with the theory of input / state / output 
linear systems for the Bergman-space setting. In particular, we show that each of 
the four approaches to Beurling-Lax representations sketched above (which blur 
into each other in the classical case) actually lead to four distinct kinds of theories 
for the Bergman-space setting. Our results for the first three approaches are new 
for the Bergman-space setting, while the fourth approach is most closely connected 
with the notion of Bergman inner function as appearing in the work of Hedenmalm, 
Olofsson, and Shimorin. Specifically, for M C A n (y), the first approach leads to 
the representation for A4 as 



j>0 



where the weights /life's arc defined by 




(1.12) 



M = M F 0^(Wj , where F{z) = [Fi{z) ■ ■ ■ F n (z 




WEIGHTED BERGMAN SPACES 



5 



and where Fj is a bounded multiplier from Aj(Uj) into A n {y) for each j = 1, . . . ,n 
with the additional property that the multiplication operator [Mp 1 ■ ■ ■ Mp n ~\ is 
a partial isometry from 0™ =1 Aj{Uj) into A n {y) (see Theorem 17. ip . The second 
approach leads to the identification of certain conditions when one has a repre- 
sentation M = MF n A n {U n ) (i.e., one can take F\, . . . , -F n -i all equal to zero in 
the previous representation) and applies to the more general situation where M. 
has its own norm and is contractively included in A n (y). The third approach 
leads to the construction of a Bergman inner family 60,61,62,- ■• of multipli- 
ers, where z k Qt{z) maps the coefficient space Uk isometrically into A n (y) so that 
[Mq Mq 1 Mq 2 ■ ■ •] maps the time-varying Hardy space 0^L O isomet- 
rically onto Ai C A n (y)- This perhaps is the most compelling new Beurling-Lax 
representation for the Bergman-space setting and has the most striking connections 
with systems theory. Namely, the multiplication operator \Mq Mq x Mq 2 ■ ■ •] 
can be identified as the Z-transform of the input-output map of a certain time- 
varying linear system (see formula Q3.2p below) having certain additional metric 
properties (see (|6.7I) , (|6.22[) and Section [Sj) . When one specializes to the classical 
case n = 1, this time- varying linear system collapses to the time-invariant linear 
system (|1.6j) with transfer function equal to the Beurling-Lax representer 6(z) in 
the classical case. 

To make these system theory connections precise, we delve into the structure of 
observability operators and observability gramians, operator resolvents, and trans- 
fer functions more general than in the classical case. Namely, observability opera- 
tors of the form 

00 

O niClA - x^C(I~ zA)- n x = J2 ( i+ J _1 ) ^ X > J ( L13 ) 

3=0 

as well as a fc-shifted version 

00 

3=0 

as well as functions having a realization of the form 

&(z) = D + z£) n ,k,c,A (z)B (1.14) 

come up. Indeed such formulas appear already in the work of Olofsson [29], at 
least for the cases k = 0, 1. We show that O nj c,A (|1-13|> arises as the observability 
operator for the time- varying system (|3.2p and that operator- valued functions of the 
form (|1.14p arise naturally in the construction of kernel functions for subspaces 
of the form S„Ai S k+1 A4 (where M. is an ^-invariant subspace of A n (y))- 

As preliminaries to the exposition of these ideas, in Section[2]we develop a calcu- 
lus of shifted power geometric series as preparation for their use in the corresponding 
operator functional calculus. In Section [3] we present the form of the time- varying 
system (|3.2p and develop its behavior under the Z-transform for the general set- 
ting where no metric constraints are imposed. Section U extends standard stability 
notions to the Bergman-setting where the system has the form (|3.2p . Section [S] 
identifies the reproducing-kernel structure on ranges of observability operators; a 
key result is Theorem 15.61 which identifies the reproducing kernel for a subspace 
of A n (y) of the form S k R&nD n ,k,c,A- The next Section |6] develops the metric 
properties for the system operators in the system (|3.2p which will be needed for the 
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Bergman-space Beurling-Lax representations to come. Finally Section [7] develops 
our Beurling-Lax representation theorems for the Bergman-space setting (all four 
approaches) while the concluding Section [5] makes the connections with the linear 
system (|3.2I) (with the additional metric constraints imposed) precise. 



2. Power-series representations for generalized geometric series 

We start with recording an assortment of power series expansions which will play 
a key role in the sequel. We let S* be the standard backward shift operator acting 
on formal power series according to 

S* : jr ajz j h> jr a j+ iz j or, equivalent^, SJ : f(z) i-> ti^hilM . (2.1) 

3=0 j=0 

Let us introduce the notation 

Rn(z) = (1 - z)- n and R n , k (z) = (S{) k {1 - Z y n , (2.2) 
so that R n ,o{z) = Rn{z), and record the power series expansions 



oo 



Rn{z) = J2 ( n4 f *) • * j > = E ( n+ ?+^~ 1 ) • ^ ( 2 - 3 ) 

The first representation follows from successive term-by-term differentiation of the 
geometric series (1 — z)^ 1 — X^>o z ^ wnereas the second follows from the first and 
the definition (|2.1j) of S{. It is seen from (|2.3p that for the special case n = 1, 

i? M (z) = i?i(z) for all A; = 0,1,... (2.4) 

The formula 

n 

R*Az) = '£ i { t+ £*)R n -i + i{z) for fe>l (2.5) 
t=\ 

follows from (12.31) and the Chu-Vandermonde identity for binomial coefficients (see 
e.g. [23 page 50]) 

n 

( n+ jir)=E(^-T 2 )-( n T')' ( 2 - 6 ) 

according to which indeed 

oo n n I oo 

*m (*) = ££( Tt 2 ) ( n T' )^' = E( ^-T 2 ) • E ( n T e ) 

j=0 £=1 t=l \j=0 

n 

where we made use of (|2.3I) again in the last step. An easy corollary of definitions 
(J2T2P is the following. 



Lemma 2.1. TTie functions R n ,k{z) given by (|2.2[) satisfy the recursion 

Rn,k{z) = ( '^r 1 ) +zi?„, fe+1 (z). (2.7) 
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Proof. Since R n M+i(z) = {Sl) k+1 (1 - z)~ n = S{ {{St) k {\ - z)~ n ) = S^Rn^z), 
we have 

R n ,k+i{z) = S{R n , k {z) = [R 7h k( z ) - Rn,k(°)]/ Z - 

From the expansion for R n k(z) in (|2.3p we see that R n .k(0) = ( n+ fc ^ an d the 
formula (12.71) follows. □ 

The following analogue of the well-known formula J2f=a z ^ = for the 

sum of the truncated geometric series will be useful in the sequel. 

Proposition 2.2. The truncation of the infinite series representation for R n (z) in 
(|2.3|) has the explicit summation formula: 

N . n N+j 

s("T 1 )=^o J i F -i:(E; , ) I fT W - < 2 -«) 

3=0 v ; 3=1 v > 

Proof. The following computation 

\ / n+j-1 » "* 1 " / ! 



3=0 



(n - 1)! dz™- 1 V 1 - 



l-z N+n 1 ^?. B _ n j!(AT + fe)! z N +J +1 



(1-2)™ J ; (iV + j + 1)! (1 - z)i +1 

(by the Leibnitz rule for the derivative of a product) 

i 7 N+k «-2 -JV+j+1 
1 z \ / N+n \ _f 

(1 - Z) k K N+j+1 > (1 - Z)3+1 

1 " 7 N+j 

verifies the result. □ 



Proposition 12.21 enables us to express the shifted n-th-power geometric series 
Rn,k(z) in terms of unshifted n-th-power geometric series R n {z) as follows. 

Proposition 2.3. The shifted n-resolvent R n ^(z) is recovered from the unshifted 
n-resolvent R n (z) according to the formula 

n-i / k \ 

R n , k (z) = R n (z)-Y, E(- 1 ) j ( fc+ r 1 )( n_1 7 K+ zK ( 2 - 9 ) 
K =o \j=o y 

Proof. By definitions (|2~2]) . 

fc-i 



i?„, fe (z) = • UnW-E("T 1 

= z- fc i?„(z) • [ 1 - (1 - z) n ■ E ( ) * j I • (2- 10 ) 



3=0 
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Letting N — k — 1 in (12. 8j) and combining this formula with (|2.10|) gives 



= ttn(Z) 


n-1 

\ " 1 k+n- 
k 

3=0 




„ \ n — J — 1 


= Rn(z) 


n — 1 n—j — 1 

E E 




("-f 1 ) (-i)V' +£ 


= Rn(z) 


n — lf k 

EE' 


; fc+ r 1 )i 




= Rn(z) 


n — lf k 

EE' 

re=0 \j=0 


; fe+ r 1 )i 


^»-i-*fi) (_i)j j z K 



and (JUnj) follows. □ 

As a consequence of identities (|2.3I) . (|2.5I) . (|2.7I) . (|2.9[) . we have the following 
identities for any Hilbert-space operator A having spectral radius less than one: 



oo 

R n {zA) = (I- zA)- n = ( B+ *~ 1 ) AV, (2.11) 

3=0 

OO 

E„ )fe (zA)=^(" + j^- 1 )^V, (2.12) 



3=0 



R n . k {zA) = Y,( i+ i k ri 2 ) R n-i+i(zA) for fc > 1, (2.13) 
1=1 

i?„, fc (zA) = (^-^Ix + zAR^+^zA), (2.14) 



n— 1 



i?„, fc M) = R n (zA) . j2 E(-!) J ( fc+ r*) (" ' ." ') (2.15) 

K=0 \j=0 J 

3. System theory motivation 

We have seen in the Introduction (see formulas (|1.6[l - (|1.9[l ) how a formula of the 
type 

6 Sl (z) = D + zC{I - zAy 1 B = D + zCR l {zA)B. (3.1) 

arises as the transfer function of a discrete-time time- invariant input/state/output 
linear system of the form (|1.6p . Furthermore, imposition of the condition that the 
colligation matrix [ q ^ ] be contractive (respectively, unitary) together with an ad- 
ditional stability condition on the state-update operator A leads to the transfer 
function being in the Schur class, i.e., having contractive values on the unit disk 
(respectively, being inner, i.e., contractive on the unit disk with unitary nontan- 
gential boundary values on the unit circle almost everywhere). 

In the quest for a parallel system-theoretic interpretation for Bergman inner 
functions on the unit disk, Olofsson [27l [29] has shown that realization formulas of 
the type Q(z) = D + zCR n {zA)B involving higher-order resolvents R n (zA) arise 
naturally and he also associated a certain higher-order time-invariant linear system 
having such a Q(z) as transfer function. We shall see below that functions of the 
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form Q(z) = D + zCR n ,k(zA)B involving shifted higher-order resolvents R n ,k{zA) 
also arise naturally. To obtain a system-theoretic connection for functions involving 
such shifted higher-order resolvents R nt k(zA), we propose to consider the following 
discrete-time time-varying linear system: 

,r ABj]] y( x(j + i) = ^ ■ Ax {j ) + ■ b M j), 
\l CD i\! j&z J ■ j y{j) = cx(j) + (^r 1 ) ■ D jU (j) 

where 

AeC(X), CeC(x,y), B k eC(U k ,x), D k e£(U k ,y) 

are given bounded linear operators acting between given Hilbert spaces X, y and 
Uk (k > 0). We note that the case where n = 1 and where the operators Bk = B 
and Dk = D are taken independent of the time parameter k G Z+ reduces to the 
classical time- invariant case given by (|1.6[) . 

If we let the system evolve on the nonnegative integers j £ Z+, then the whole 
trajectory {u(j),x(j),y(j)}j^z + is determined from the input signal 
and the initial state x(0) according to the formulas 



<j) - ( "T 1 ) ■ A3 <°) + E A-^Bmie) , (3.3) 
= ( "+f 1 ) ■ (cA>x(0) + £ CAi-t-'BMt) + D 3 u(j)) ■ (3-4) 



1=0 

Formula (|3.3[) is established by simple induction arguments, while (I3.4[) is obtained 
by substituting (|3.3|) into the second equation in (|3.2|l . 

To write the Z-transformed version of the system-trajectory formula (I3.3p . we 
multiply both sides of (|3.3p by z J and sum over j > to get, on account of (|2.11l) . 

oo 

x(z) = ^x{i)z 

3=0 

D \ OO / OO \ 

(^r 1 ) AV *( Q ) + E E ('T 1 ) - !) 

v i=o y fe=i y=/s y 

CO / oo \ 

(/ - z^-^o) + E ^ E ) AV - x ) 
fc=i \j=o y 

oo 

(/ - zA)- B a:(0) + ^ z k R n , k {zA)B k -iu(k - 1) 
fe=i 

oo 

(/ - zA)- n x(0) + z k+1 R n ,k+i(zA)B k u(k). (3.5) 



k=0 
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The same procedure applied to (|3.4|) gives 

oo 
k=0 

oo 

= C(I - zA)- n x(0) + J2 zk (( fe+ r' ) D k + zCR n , k+1 (zA)B k ) u(k) 

k=0 

oo 

= O n ,c, A x{0) + **©»,*(*)«(*), (3.6) 

fc=0 

where 

oo 

O n , C ,A: x^Y, (("T 1 ) CA3x ) z ' = C ( J - xA )~ nx ( 37 ) 

is the n-observability operator and where 

e n ,k(z)=( k +l- 1 )D k +zCR n , k+1 (zA)B k {k = 0,1,...) (3.8) 

is the family of transfer functions. We note that the transfer function Q n ,k(z) 
encodes the result of a pulse input- vector u being applied at time j = k: 

y(z) = ® n}k (z) ■ z k u if x(0) — and u(j) = 8^ k u. 

In fact the functions Q n ,k(z) could have been derived in this way and then one could 
arrive at input-output relation (|3.6I) via superposition of all these time-/c impulse 
responses. Note also that formula (|3.8p for the classical time-invariant case (n = 1 
and B k , D k independent of k) reduces to the formula (|3.1|) for the classical transfer 
function Qs 1 due to the identity (|2.4j) . 

For the application to Bergman inner functions, it is natural to impose some 
additional metric constraints on the colligation matrices [ q p k ] ; these are discussed 
in Section [HI below. 

4. Observability operators and gramians, Stein equalities and 

inequalities 

Formula (|3.7| associates with any output pair (C, A) (i.e., C £ C(X,y) and 
A e C{X)) the n-observability operator O n ,c,A- In case O n ,c,A is bounded as an 
operator from X into An(y), we say that the pair (C, A) is n-output-stable. If 
(C,A) is n-output stable, then the n-observability gramian 

Gn,C,A '■= (Cn : C : A)*Cn,C,A (4.1) 

is bounded on X and can be represented via the series 

oo 

Gn,c,A = ( fe+ r' ) A* k C*CA k (4.2) 

k=0 

converging in the strong operator topology. As suggested by the Agler hereditary 
functional calculus as formulated by Ambrozie-Englis-Miiller [4], we introduce the 
operator 

B A :X^A*XA (4.3) 
mapping C(X) into itself, and then view GnC A (at least formally) as being given 

by 

Gn,c,A = (I-B A )- n [C*C]. (4.4) 
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If \\Ba\\ < 1, (|4.4|) is precise; in general one can make this precise by interpreting 
(|4.4j) in the form 

g n ,c,A = lim(/ - rS A )-"[C*C]. (4.5) 

rfl 

Either of the formulas (|4.4[) and (|4. 5[) suggests that we define 

Go,c,a = C*C. (4.6) 
We next introduce the operator map 

T n , A = (I-B A ) n : X^^{-l) k ( n \A* k XA k . (4.7) 

k=0 ^ ' 

There is then an assortment of identities as listed below. 

Lemma 4.1. (1) For all H,A^ C(X) and any integers k > 1 and N > 0, 
T k ,A[H] = r k - ltA [S] - AT fc _i, A [fl]A, (4.8) 

N k 

H = J2 ( k+ i~ ') A*lT kiA [H]Ai ! + (%Xj) A* N ^T^ jtA [H}A N ^. (4.9) 

3=0 j=l 

(2) If (C, A) zs n-output stable, then also 

Q m ,c,A- A*g m! c,AA = g m -i t c,A for m = l,...,n, (4-10) 

r^A^n.C.A] = g n -k,C,A for k = 0,...,n. (4-11) 

Moreover in this case (C, A) is also k-output stable for k = 1, . . . , n — 1. 

Proo/. Upon applying identity (J - = (7- B A ) k ^ - B A (I - Ba)^ 1 to an 

operator H e C(X) and making use of definition (|4.7p we get (|4.8j) . To verify (|4.9|) . 
we use the identity (|2.8|) in the form 

N JV+7 

\ ^ / fe+i-i A „j _ 1 _ ( N+k\ J_ 

3 ) - (i _ z) k n+3 ) (1 _ z y ■ 

3=0 y ' j=l V > 

Multiplying both parts in the latter identity by (1 — z) k we get 

N k 

1 = E ( ^ ) ^ - *)* + E ( nX" ) *" +J (l " z) N ~ j i ( 4 - 12 ) 



3=0 3 = 1 

which in turn implies the operator identity 

JV A 
3=0 j=l 

Upon applying this latter identity to an operator H 6 £{X) and making use of 
definition gjj) we get (|4T9j) . 

We provide two proofs of (|4.10[) and (|4.11l) . one of which is a straightforward 
easy-to-remember computation but which is rigorous only if the spectral radius 
p(Ba) of B A is less than 1, together with a second more elaborate proof to handle 
the general case. We have from f|4.3|) and (|4.4[) 

Gn.C.A - A*g m , C ,AA — {I - BA)[Gm,C,A] 

= (I-B A )o(I-BAT m {C*C} 

= (i - B A y m+1 [c*c] = g^.c.A 
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and (|4.10p follows, at least for the case where \\Ba\\ < 1. To handle the general case 
(where we only assume that (C, A) is n-output stable) , one can plug in the infinite 
series representation (|4.2p for Q m ,c,A and make use of the binomial coefficient iden- 
tity (™) = ('"jT 1 ) + {"kZx) to arrive at the result. A corollary of the computation is 
that the infinite series defining G m -i,C,A is strongly convergent, i.e., (C, A) is also 
(to — I )-output stable whenever it is m-output stable. 

To prove (|4. 1 If) under the assumption that p{Ba) < 1 , note that 

Tj, A [Gn,c,A] = {!- b a) j [(I - B A )- n [C*C]] 

= {I- B A )- n+ i[C*C] = g n - j>c ,A- 

To handle the general case, one can do an inductive argument using identity (|4.f 0[) 
to arrive at the result. As a corollary we arrive at the final statement in Lemma 
14. It (C, A) is fc-output stable for k — l,...,n — 1 whenever (C, A) is n-output 
stable. □ 

Definition 4.2. The operator A e C(X) is called n-contractive iiT n}A [I] is positive 
semidefmite and it is called n-hypercontractive if T^ _a[J] > for all < k < n. 

It was shown in [1] (see also [25] as well as [26] for a multivariable version) that 
inequalities r^f/] > and r„ ^[7] > imply that A is an n-hypercontraction. 
This result extends from / to an arbitrary H > 0; the proof below is modelled from 
the one in [25] . 

Lemma 4.3. Let us assume that the operators H, A € C(X) are such that 

H>A*HA>0 and T n>A [H] > (4.13) 
for some integer n > 3. Then 

r k , A [ H }>° for all k = l,...,n-l. (4.14) 



Proof. Observe that the leftmost inequalities in (|4.I3I) mean that H and Ti^H] 
are both positive semidefmite. Making use of definitions (14.31) and (14. 7|) we have 

N N 

Y,A* j r n _ hA [H]Ai =Y, B aV 

3=0 3=0 

= (I - B% +1 )(I - B A ) n - 2 [H] 

= T n ^ A [H] - A* N+1 T n ^ A [H]A N+1 . (4.15) 
Iterating the first inequality in (|4.I3|) gives A* 1 HA 3 < H for all j > and therefore, 

n-2 

\(T n ^ A [H]A k x, A k x)\ < (V) ( Hx > x ) = 2 "~ 2 <^ *)■ ( 4 -!6) 

3=0 

Taking the inner product of both parts in (14. 15|) against x € X and then making 
use of (|4.16p gives 



A' 



J2 {^n-i AH}A 3 x, A°x) 

3=0 



= \(T n ^ A [H]x, x) - (T n ^ A [H}A N+1 x, A N+1 x)\ 
< 2"- 1 (Hx, x) . (4.17) 
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On account of relation (|4.8j) . the second inequality in ()4.13|) implies 

T n -i, A [H} > A*T n _ ltA [H]A. 

Therefore (T n _ ltA [H]A>x, A>x) > (T n ^ A [H}A^ +1 x, A> +1 x) for all j > and 
all a; S X. Thus, on the left hand side of (14. 17)) we have the partial sum of a 
non-increasing sequence and, since the partial sums are uniformly bounded (by 
2«-i (ffx, x)), it follows that all the terms in the sequence are nonnegative. In 
particular, (T n -i >A [H]x, x) > 0. Since the latter inequality holds for every x G 
X, we conclude that T n -i tA [H] > 0. We then obtain recursively all the desired 
inequalities in (|4.14|) . □ 

Definition 4.4. The operator A 6 £>{X) is called strongly stable if A k tends to 
zero as k — > oo in the strong operator topology, i.e., ||^4 fe a;|| — > for every x 6 X. 

The following result gives connections between n-output stability, observability 
gramians and solutions of associated Stein equations and inequalities. In case n = 1 
it is well-known. In what follows, we will refer to the last relations in (|4. 18[) and 
(|4.19l) as the Stein inequality and the Stein equality, respectively. 

Theorem 4.5. Let C G £(X,y) and A G C(X). Then: 

(1) The pair (C, A) is n- output- stable if and only if there exists an H G £{X) 
satisfying the system of inequalities 

H > A* HA > and T„ jA [H] > C*C. (4.18) 

(2) // (C, A) is n- output- stable, then the observability gramian G n ,C,A satisfies 
the system 

H > A* HA > and T n>A [H] = C*C (4.19) 

and is the minimal positive semidefinite solution of the system (|4. 18[) . 

(3) There is a unique positive- semidefinite solution H of the system (|4.19|) with 
H = Q n ,c,A if A is strongly stable. If A is a contraction, then the solution 
of the system (|4.19|) is unique if and only if A is strongly stable. 

Proof. Suppose first that (C, A) is n-output-stable. Then the infinite series in (14.2[) 
converges in the strong operator topology to the operator H = Q n .c. A > and by 
(1440)1 . f4~TTT) and (l4H . 

Qn,C,A — A*Q n ,c, A A — G n -1,C,A > 0, T n ^ A [Qu.c.a] = Qo,C,A = C*C. 

Thus, H = Qn.c.A satisfies relations (14.19[) and hence also inequalities (|4.18l) . 

Conversely, suppose that inequalities (|4. 18[) are satisfied for some H £ C(X). 
Then for every integer N > we have 

N N 

J2 (^+Y 1 ^A* 1 C*CA 1 < A* j T n , A [H]A> 

3=0 j=0 

n 

= H- (%X1) A * N+Jr n-j, A [H]A N+ i < H. (4.20) 

3=1 

Indeed, the first inequality follows since C*C < T n , A [H], the second equality holds 
by (|4.9|) (with k = n) and the last inequality holds since r^fiS] > for k = 
0, . . . , n — 1, by the assumptions (|4.18|l and Lemma [43l 
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By letting N — > oo in (|4.20[) we conclude that the left-hand side sum converges 
(weakly and therefore, since all the terms are positive semidefinite, strongly) to a 
bounded positive semidefinite operator. By (|4.2[) . 

N oo 

lim A * jC * CAj = J2 ( i4 7 _1 ) A * iC * CAi = Gn,C,A 

3=0 ' 3=0 

and passing to the limit in (|4.20|) as N — >• oo gives £7 n ,c,A < -ff- In particular the 
operator 5n,c,A is bounded (since H is) and therefore the pair (C, A) is n-output- 
stable. This completes the proof of the first part of the theorem. 

As observed above, the gramian Gn.c.A satisfies relations (|4.19|) (and therefore, 
relations (|4.18|) ). We also showed that any operator H solving the inequalities 
f|4.18|) also satisfies the inequality H > Q n ,c,A so that Q n ,c,A is indeed the minimal 
solution to the system of inequalities (14. 18)) . This completes the proof of the second 
part of the theorem. 

Now suppose that A is strongly stable and that H G C{X) solves the system 
(|4.19[) . We will show that then necessarily H = Q n ,C,A- We first recall that if 
positive semidefinite operators P, Q £ C(X) satisfy the Stein equation 

P-A*PA = Q (4.21) 

with a strongly stable A G C(X), then P is uniquely recovered from (|4.21[) via the 
strongly converging series 



P =y^A* k QA k . (4.22) 

Indeed, iterating (|4.21j) gives 



k=0 



N-l 

P = A* N PA N + A * j Q Aj for all N > 0, 

3=0 

and since all the terms in the latter equality are positive semidefinite, the conver- 
gence of the series on the right side of (|4.22l) follows. The strong stability of A 
guarantees that A* N PA N ->• as N ->• oo, which implies equality in (I4.22[) . Now 
we observe that by (|4.8[) . the Stein equation in (|4. 19[) can be written as 

r n _ 1>A [#] - A*T n _ ltA [H]A = C*C 

and the above uniqueness shows that 

oo 

r„_ M [H] = A * k c*CA k = g h c, A . 

The latter equality can be in turn written in the form (I4.21[) with P = T n _2A[H] 
and Q = Q\ t c,Ai an d by the same uniqueness argument we have 

oo 

r n -2,A[H] = ^ A* k Qi l c,AA k — G2,C,A- 

Continuing this procedure, we get r n _ [iJ] = Gj.c.A for j = 1, . . . ,n. For j = n 
we have in particular, H = Tq,a[H] = Q n ,c,A, which gives the desired uniqueness. 

To complete the proof of part (3) of the theorem it remains to show that if A is 
a contraction and the system (14. 19)) admits a unique solution (which necessarily is 
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H = Qc.A.n), then the operator A is strongly stable. We prove the contrapositive: 
if A is not strongly stable, then the solution of (|4.19|) is not unique. 

Since A is a contraction, the sequence of operators A& = A* k A k is decreasing 
and bounded below and therefore has a strong limit A > which clearly satisfies 
the relation Ba[A] = A and which is not zero, since A is assumed not to be strongly 
stable. Then it follows that 

k 

T k<A [A] = (I- B A ) k [A] = (*) A = 

j=o 

for all k G {1, . . . ,n}. Now we see that the operator H — Gn.c.A + A (as well 
as Gu,c,a) satisfies the system (|4.19[) which therefore has more than one positive- 
semidefinite solution. □ 

Definition 4.6. A pair (C, A) with C € C(X,y) and A £ C(X) is called n- 
contractive if inequalities (|4.18l) hold with H = Ix, that is, if A is n-hypercontractive 
and in addition 

n 

T„,a[4] > C*C, or equivalents C*C + £ (" ) < 4 . 

i=i 

The pair (C, A) will be called n-isometric if relations (|4.19l) hold with H = Ix, that 
is, if A is n-hypercontractive and in addition 

n 

r„, A [Ix] = C* C, or equivalents C* C* + ^ ( ™ ) A* J A J = 7 A , . 

i=i 

A pair (C, A) is called observable if the operator Ox,c,A (equivalently, ^1,0,^1) 
is injective. This property means that a state space vector 1 £ ^ is uniquely 
recovered from the output string {yk}k>o generated by running the system (13. 2[) 
with the initial condition xo = x and the zero input string. A pair (C, A) is called 
exactly observable if Oi t c,A (equivalently, Gi,c,a) is bounded and bounded from 
below. More generally, we say that the pair (C, A) is exactly n-observable if the 
rt-gramian Q n ,c,A is bounded and bounded below. While n-observability implies k- 
observability for 1 < k < n, the corresponding statement for exact fc-observability 
fails (see Proposition 15. 71 below). 

The following statements follow along the lines of Theorem 14.51 

Proposition 4.7. (1) Suppose that (C, A) is an n-contractive pair. Then (C, A) is 
n- output- stable with Q n ,c,A < I x and the observability gramian Q n ,c,A is the unique 
solution of system (|4.19p if and only if A is strongly stable. 

(2) Suppose that (C, A) is an n-isometric pair. Then H = I x is the unique 
solution of the system (|4.19j) if and only if A is strongly stable. In this case O n ^c,A 
is isometric and hence also (C, A) is exactly observable. 

Let us say that the pair (C, A) is similar to the pair (C , A) if there is an invertible 
operator T 011 X so that C = CT- 1 and A = TAT' 1 . Then we have the following 
characterization of pairs (C, A) which are similar to an n-contractive or to an n- 
isometric pair. 
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Proposition 4.8. (1) The pair (C, A) is similar to an n- contractive pair (C, A) if 
and only if there exists a bounded, strictly positive-definite solution H to the system 
of inequalities (|4.18l) . 

(2) The pair (C, A) is similar to an n-isometric pair if and only if there exists 
a bounded, strictly positive-definite solution H of the system (|4.19p . 

Proof. Suppose that H is a strictly positive-definite solution of (|4.f 8|) . Factor H 
as H = T*T with T invertible and set 

C=CT~ X , A^TAT- 1 . (4.23) 

Multiplying each inequality in (|4.18j) on the left by T*^ 1 and on the right by T^ 1 
and replacing A and C respectively by T~ X AT and CT then leads us to 

I X >A*A and T nA [I x ]>C*C 

which means that (C, A) is an n-contractive pair which is similar to the original 
pair (C, A). Conversely, if (C,A) given by (|4.23j) is contractive, then H = T*T 
is bounded and strictly positive-definite and satisfies the inequalities l|4.18jl . This 
verifies the first statement of the Proposition. The second statement follows in a 
similar way. □ 

As a consequence of the observations in Proposition 14. 8[ Proposition [477] can be 
formulated more generally as follows. 

Proposition 4.9. (1) If the pair (C, A) is such that inequalities (|4.18j) have a 
strictly positive-definite solution H, then (C, A) is n- output- stable. Moreover, the 
observability gramian G n ,c,A is the unique positive semidefinite solution of the sys- 
tem (|4.19[) if and only if A is strongly stable. 

(2) If the pair (C, A) is such that the system (|4. 19|) has a strictly positive-definite 
solution H , then (C, A) is n- output- stable and the observability gramian Q n ,c,A is 
the unique positive semidefinite solution of the system (14. 19)) if and only if A is 
strongly stable. In this case (C, A) is moreover exactly observable. 



The last part of Proposition 14.91 has a converse. 

Proposition 4.10. Suppose that the pair (C, A) is n- output- stable and exactly 
observable. Then A is strongly stable. 

Proof. Observe that if H is any solution to the system (|4.19|) , then (|4.20|) takes the 
form 

N n 

) A * jc * CAj = H -J2( nVj ) A* N+j T n _ jA [H]A N+j (4.24) 

J=0 3=1 

and still holds for all N > 0. If (C, A) is n-output-stable and exactly observable, 
then the observability gramian Gn.c.A is a strictly positive-definite solution of the 
system l|4.19[) . Plugging H = Q n ,c,A into (I4.24|) gives 

N n 

G n ,c,A =J2 C* 4 *" 1 ) A* k C*CA k + J2 A* N+ lr n _ jjA [g n ,c, A }A N+ i 
k=0 j=i 

N n 

= Y,C*V 1 ) A "' C * CAk + E (% + +1) A * N+3 Sj.c,aA N+3 (4-25) 
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where the last equality follows from (|4. 1 lj) . From the infinite-series representation 
(|4.2|) for Gc,A,n, taking limits in (|4.25|l gives 

Um (»+n} A *N +j g A N +j = Q 

for j = 1, . . . , n. In particular, we have 

lim A* N+n g n ,cM N+ " = 0- (4.26) 
Since G n ,C,A is strictly positive definite, we conclude that there is an e > so that 

£\\x\\ 2 < {Gn,c,A%, x) for all x e X. 
Upon combining (|4.26|) with the latter relations (with x replaced by A N x) we get 

s\\ A N xf < (G n ,cM N x, A N x) 
for all x £ X, and we conclude that A is strongly stable as asserted. □ 

We conclude this section with formal introduction of backward-shifted versions 
of the observability operator O ni c,A (|3.7I) and observability gramian G n ,C,A (|4.1[l . 
The latter operators can be expressed as 

O n<c ,A ■ x -> CR n (zA)x and Q nfi ,A = R n (B A )[C*C] 

where R n (z) = (1 — z)~ n and where Ba is the operator on C(X) given by (I4.3p . We 
now introduce the backward-shifted variants of these objects by simply replacing 
the function R n in the two formulas above by its backward shifts R n ,k- 

oo 

O n ,k,c,A : x -> CR n , k (zA)x = ( n+ ££ _1 ) {CA^x)z\ (4.27) 

3=0 

OC 

®n,k,c,A := JZ»,*(Ba)[C*C7] = E ') ^*C*CA'. (4.28) 

Observe that the latter power series representations follow from that in (|2.12l) . 
Letting k = in flj^TTJ) , (|4~2"gf we conclude 

O«,0,C,A = Cn,c,A and ©„,o,C,,4 = G n ,C,A- (4.29) 

If fc > 1, then we can use the formula (12.5[) along with (|4.27|) . (I4.28|) to get repre- 
sentations 

n n 

&n,k,C,A = 2^ ( £ ^-l 2 ) ^n-<+l,C,j4) ©Ti,fc,C,A = 2_, ( £ ^e— T 2 ) &»-*+l,C,.A 

e=i t=i 

which imply in particular that for an n-output stable pair (C, A), the operator 
n,k,c,A ■ X — > A n (y) is bounded for all k > 0. 

Proposition 4.11. The weighted Stein identity 

A*® n<k+l< cM + ( n+ fe _1 ) • C*C7 = <8 n , fe , c , A (4.30) 
ZioZds /or all integers k > 0. 
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Proof. The operatorial equality 

R n ,k(B A ) - B A o R n .k+i{B A ) — R n ,k(0) ■ I c(x) — ( n+ fc _1 ) ' Ic(x) 

is a consequence of the identity (|2.7I) . Applying this equality to the operator C*C 
gives, on account of (|4.28[) . 

(R n ,k{B A ) - B A o R n . k+1 (B A )) [C*C] = & n , k , c ,A ~ A*<8 n , k+hC ,AA 

= ( n+ ^ 1 ) ■ C*C 

and we arrive at (|4.30[) as wanted. □ 

As a consequence of Proposition |2~31 it turns out that & n k,c,A can De expressed 
in terms of Q n ,c,A a $ follows. 

Proposition 4.12. The operator <& n ,k,c,A given by ()4.28j) admits the representa- 
tion 

n—l/K \ 
®n,k,C,A = E E ( ) ( ) ^"Sn.C.AA". (4.31) 

Proof. Note that the identity (12. 9[) leads to the operator identity 
n—X/K \ 

w = E E(T) (" '/ ') (^roi?„(s A ). 

K=0 \j=0 J 

Application of this operator to C*C now leads to the identity (|4.31|) due to (|4.28[) . 
(IQll and gH. □ 

5. Observability-operator range spaces and reproducing kernel 

HlLBERT SPACES 

Let S n denote the shift operator on the space A n (y) defined as S n : f(z) — > 
zf(z). Reproducing kernel calculations show that its adjoint S 1 * is given by 



Iterating the latter formula gives 

o*fcf V' (k+j)\(n + j - 1)! ^ Mnj+fc f j ,, 9 n 

5 « / = g ,!(n + ^-l)! ' ^ - g ^7 ' ■ ^ 

The operator S 1 * is a strongly stable n-hypercontraction. As was shown in [28l 
Lemma 5.1] for / of the form (|5.1I) . 



E(-i) J ll^/ll; 



= ll/of ■ (5.3) 



The latter equality means that the pair (E.S*) is n-isometric in the sense of Defi- 
nition 223 where E : A n (y) — 5- 3^ is the evaluation operator defined by Ef = /(0). 
It follows from flOl) that 

nHn-lV 

ES n f = /lu = ^,kh (5.4) 

(n + k ~ 1)! 
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and therefore, 

oo oo 

O n ,E,Sif = E(I - zS* n )- n f = ( J+ T 1 ) ( ES nf) * = E = f 

3=0 j=0 

so that the observability operator O n ^E,s* equals the identity operator on A n (y)- 
Associated with an n-output-stable pair (C, A) is the range of the observability 
operator 

RimO n ,c,A = {C(I - zA)~ n x: x e X}. 

Theorem 5.1. Suppose that (C, A) is an n- output- stable pair. Then: 

(1) The intertwining relation 

S* n O n>c ,A = O n , c , A A (5.5) 

holds and hence the linear manifold M = Ran0 n ,cA is S* -invariant. 

(2) Let H £ C(X) satisfy inequalities (|4.18[) and let X 1 be the completion of 
X with inner product \\[x\\\ 2 x , — {Hx,x)x (where [x] denotes the equiva- 
lence class modulo Ker H generated by x). Then A and C extend to define 
bounded operators A' : X' — > X' and C : X' — > y and the observability 
operator O n: c,A extends to define a contraction operator O n ,c'.A' '■ X' — > 
A n (y). Moreover, O ni c',A' * s an isometry if and only if H satisfies rela- 
tions (|4.19[) and A' is strongly stable, i.e., (HA N x, A N x) — > for all x £ X 
as N — >• oo. 

(3) If the linear manifold A4 := Ka,nO n> c,A is given the lifted norm 

\\O n , c ,Ax\\ 2 M = inf {Hy,y)x, 

y€X: O n:C ,Ay=O n ,c.AX 

then 

(a) A4 can be completed to M! = RanC^c'.A' with contractive inclusion 
in A n {y): 

ll/lll < \\f\?M> f° ral1 fzM'. 
Furthermore, M.' is isometrically equal to the reproducing kernel Hilbert 
space with reproducing kernel K n ^c,A,H given by 

K n , c ,A, H (z, C) = C(I - zA)- n H(I - (A*r n C*. (5.6) 

(b) The operator S* : M. — > Ai defined as in ()5.ip is a contraction and 



(j) \\S* j ff M > 11/(0)115- for all feM. (5.7) 



y^_iY7' I n I IIO*7fl|2 || f(n\\\2 
3=0 

Moreover, if the Stein equality in (|4.19[) holds, then (|5.7[) holds with 
equality. 

(4) Conversely, if M is a Hilbert space contractively included in A n (y) which 
is invariant under S* which in turn is a contraction on M. and for which 
(|5.7|) holds, then there is an n-contractive pair (C, A) such that M = 
Ti(K n ,c,A,i) — Ran Bi c,i isometrically. In case (|5.7|) holds with equal- 
ity, then (C, A) can be taken to be n-isometric. The canonical-model choice 
of such a pair (C, A) is 

(C,A) = (E\ M ,S*\ M ) where E: f £ A n (y) ^ /(0). 
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Proof of (1): Making use of power series expansion (|2.3|) and of (|5.1|) wc get (|5.5[) : 

S:O c ,A, n x = S* n (C(I - zA)- n x) = f; £±I • h + J) (CA^ 1 *)** 

fc=0 n ^ ' 

= C{1 - zA)- n Ax = O c ,A, n Ax. 

Proof of (2): Inequalities (I4.18[) amount to the statement that the pair (C, A) 
is n-contractive and well-defined on the dense subset [X] of X' (where [x] is the 
equivalence class containing x) and hence extends to an n-contractive pair (C, A 1 ) 
on all of X 1 and moreover relations (|4.20l) hold for all N > 0. By Theorem 14.51 
G n ,c,AX < H and therefore 

\\On,C,Ax\\A n (y) = (Gn,G,AX, x) X < (Hx, x) X , 

so O n ,c,A is contractive from X to A n {y). It follows from (I4.20[) that the limit 



A H , A = Vw Y, ( N+j ) A* N+ iT n _ jA [H]A N +i (5.8) 
j=i 

exists in the strong sense. We have from (|4.20l) 

N 

< ]T ( i+y 1 ) A* j (T n , A {H] - C*C) A? 



3=0 

n N 

-~ H -J2( nVj ) A' N+ Sr n - jt A [H]A N +i - ) A * j C*CA j 

J = l 3=0 

and letting N — > oo we get, on account of (|5.8[) and 



< A * j ( r «,A[-ff] - C*C) A> <H- A H , A - g ntC ,A. (5.9) 

3=0 

By definition, G n ,C>,A>'- X' — > A n {y) being an isometry means that Q n ,C,A = H 
in which case (|5.9p forces Ah, a = and equalities throughout (|5.9|) . Since all 
the terms on the right hand side of (|5.8[) are positive semidefinite, the condition 
A//^ = implies in particular that 

0= lim A* N+n T A [H]A N+n = lim A* N+n HA N+n = lim A* N HA N , 

N— >oo ' N—>oo N— >oo 

so that A' is strongly stable. All the terms in the series in (|5.9p are nonnegative 
and therefore, each term equals zero. The term corresponding to the index j = is 
T n , A [H] - C*C. Hence T n , A [H] = C*C and thus the Stein equality in (|4~19| holds. 
Conversely, if A' is strongly stable and relations (|4.19j) hold, we conclude as in the 
proof of Theorem 14.51 that Q n ,C,A' = H, i.e., that O n ,c',A' is an isometry from X' 
into A n {y). 

Proof of (3a): Statement (3a) follows from general principles laid out in [5] (see also 
[TUHH] for applications very close to the context here). For the sake of completeness 
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and since the same construction arises again in the sequel, we sketch the argument 
here. Note that for x G X and y (z y we have 

((O n<c< Ax) (0,y)y = (C(I - (A)- n x,y) y 

= (x,H(I-(A*y n C*y) x , 

= (C(I - -A)- n x, C(I - -A)- n H(I - (A*)- n C*y) M 

where the last step follows from the assumption that O n ^c,A is isometric from X' 
into M. . From this computation we see that 

((On.c.Ax) (C),V)y = {On.C.AX, K n . C ,A,H{- ,QV) M 

and it follows that K n c.a.h is the reproducing kernel for M. as asserted. 
Proof of (3b): For / of the form f(z) = C(I - zA)~ n x, we have 

Now it follows from (15.51) that 



\\S?f\\HlK»,a,*,B) = ( HA3x > A ' x )x fOT j^ 1 - 
With these substitutions, we see that S* is a contraction on M. if and only if 
(Hx,x)x > (HAx, Ax)x for ah % G X or, in operator form, H > A* HA, and that 
on another hand, inequality (|5.7p is equivalent to 



J2(-iy ( k ) {HA 3 x, A j x) x > \\Cx\\ y for all x G X, 



or in operator form, 



j>2(~iy ( k ^j A* j HA 3 > C*C (5.10) 



3=0 v ' 

with equality in (|5.7[) equivalent to equality in (|5.10l) . To complete verification of 
part (3b), it remains to remark that relation (j5.10[) coincide with that in (|4.18p (or 
with that in (|4.19p in case equality holds in (I5.10[l ). 

Proof of (4)'- Suppose that M. is a Hilbert space included in A n (y) which is 
invariant under S* which is contractive on j\4 and let us assume that the in- 
equality ([517]) holds. Let A = S* n \ M and let C be defined by Cf = /(0) for 
all / G M. In other words, C = E\ M - Inequalities S n S* < Im and (15. 7j) 
mean that the pair (C, A) defined this way is n-contractive. Moreover, if equal- 
ity holds in (|5.7[) for every / G M, then the pair (C, A) is n-isometric. It is 
readily seen that O n .c,A = O n ,E,s^ \m = Im- Therefore, for each / G M. we have 
\\f\\u(K n c a i) — WJWm an d thus M — T-L(K n ^c,Aj) isometrically. It then follows 
from part (3a) of the theorem that in fact A4 is contractively included in A n (y) ■ O 

The following corollary is a simple, useful special case of the general situation 
laid out in Theorem 15.11 

Corollary 5.2. Suppose that the n-output stable pair (C, A) is also exactly n- 
observable. Then there is an S* -invariant subspace A4 C A n {y) so that (C, A) is 
similar to the n-isometric pair (C, A) given by 

{C,A) = {E\ M ,S* n \ M ), where E: f G A n {y) ^ /(0). 
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Proof. One can simply check directly that O n ,c,A- X — > A n {y) considered as a 
linear transformation from X onto M := KaiiOn,c,A provide the required similarity 
transformation. □ 

As explained by part (4) of Theorem 15. 11 for purposes of study of contractively- 
included, S* -invariant subspaces of A n (J 7 ) which satisfy inequality (|5.7|) , without 
loss of generality we may suppose at the start that we are working with X' as the 
original state space X and with the solution H of inequalities (|4.18p to be normalized 
to H = lx- Then certain simplifications occur in parts (l)-(4) of Theorem 15.11 as 
explained in the next result. 

Theorem 5.3. Let (C, A) be an n- contractive pair with C G C(X,y) and A G 
C(X). Then: 

(1) (C, A) is n- output- stable and the intertwining relation (|5.5p holds. Hence 
Ran O^c^ is S* -invariant. 

(2) The operator O n , c ,A ■ X ->■ A n (y) is a contraction. Moreover On.c.A is 
isometric if and only if (C, A) is an n-isometric pair and A is strongly 
stable. 

(3) // the linear manifold M := R&nO nt c,A is given the lifted norm 

\\O n ,c,Ax\\ M = HQzll* (5-11) 

where Q is the orthogonal projection of X onto (Ker O ny c,A) ± , then H} c,A 
is a coisometry of X onto A4 . Moreover, M. is contained contractively in 
A n (y) and is isometrically equal to the reproducing kernel Hilbert space 
T~L{Kn t c,A) with reproducing kernel K n _c,A{z, C) given by 

K n ,cA*, = Ci 1 - zA)~ n {I - (A*)- n C*. 

(4) If M — Ran (D^ca = H(K„,c,a) is given the lifted norm || • \\-H(K n c a) as 
in (|5.1ip . then S^\m is a contraction and 

n 

E^ 1 )' (*) ■ W S n j f\\ 2 n(K n ,a, A ) > If Will for all f G H(K n , c , A ). (5.12) 

Moreover, (|5.12p holds with equality if and only the orthogonal projection 
Q of X onto (Ker O ti ^c,a) 1 ' is subject to relations 

Q > A*QA and T nA [Q] = C*C. (5.13) 

In particular, if (C, A) is observable, then (|5 . 12[) holds with equality if and 
only if (C, A) is an n-isometric pair. 

Proof. Statements (l)-(3) and all but the last part of statement (4) are direct 
specializations to the case H = Ix of the corresponding results in Theorem 15.11 It 
remains only to analyze the conditions for equality in (|5 . 12|) . 

From the intertwining relation (|5.5p . we see that inequality S n S* < Im and the 
equality in f)5 . 12[) for a generic element / = O n ,c.AX G A4 mean that 

\\^n,C,AX\\u(K niCiA ) > \\O n ,C,AAx\\ H { Kn C A ), 
n 

]T(-ir ($) • \\O n ,c,A^x\\ 2 niKn ^ A) = \\Cz\\* 

3=0 
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for all x S X. By the definition (|5.1ip of the %(2f n> c,A)-norm, the latter relations 
can be written as 

n 

\\Qx\\ x > \\QAx\\x and '52(-l) S (j)-\\QA i z\\x = \\Cx\\l 

3=0 

which in turn is finally equivalent to (|5.13l) . □ 

Remark 5.4 (Variations on a theme: shift-invariance of the range of 
an observability operator). A fundamental observation for us is the fact that 
KanO n ,c,A is invariant under the 5* due to the intertwining condition (j5.5[) . How- 
ever, for the case of a backward-shifted observability operator O n> k,C,A, we see that 
HsLTi£) n k c,A is not S* -invariant; a partial substitute is the relation 

Sl&n,k,C,A — &n,k+l,C,A-A 

which follows directly from (12.141) . We now list some additional curious connections 
between the ranges of S*&n,k,C,A an d the operators S n and £*; we shall not have 
need of any these in the sequel. 

Proposition 5.5. If the pair (C, A) is n-output stable, then the following relations 
hold: 

SZS*D n , k , c ,A=S*- l On,*-i,c,A for all k>l. (5.14) 

Q*mcki~\ / Sn &n,k— m,C,A if m < k, ■. 

!>n * n Vn,k,C,A - j O n ,c, A A m - k if m > k. [bAb) 

S„(5*5 n )- 1 ^- 1 £)„, fc _i ) c,A = S*D n , k , c ,A for all k>\. (5.16) 
Proof. To prove ([5~T4]) . we note that by ([23|) . 



OO 

(S*i2»,*)(*) = £ ( 



n+j+fc-l ^ _ z j + k 



j+k 

3=0 



and therefore, we have from (15.1 



S* S k B u(z) - V k+J fn + k + j - l\ k+ 



3=0 

OO y s 

= £ r t+ 1 - 1 ) • * fc+J ' _1 = ^^-ito - ^i^fc-iW 

which together with (|4.27|) implies (15.141) . 

The identity (j5 . 15[) follows recursively from (15. 14|) and (|5.5p . Note that the 
special case k = 0,m = 1 gives us back (|5.5I) . To verify (j5.16[) . we first note that 



J=0 j=0 ' 

from which it follows that S*S n is invertible on ^4„(3^) and moreover, 



J=0 j=0 



Mss^r 1 : E>' • -Ej-r 
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We now note that (|5T6|) is an immediate consequence of (|5.14j) : 

'A 

5 n D„ ~ZJn,k,C,A = d n *->n,k,C,A- 



S n (S^S n ) 1 &n,k-l,C,A — S n {S^S n ) SnS n £> n ,k,C,J 

= SnSt^On.k.C.A = S k O r 



□ 



Finally we note that the operator S n (S*S n ) was introduced and studied in 
the general context of a left-invertiblc Hilbcrt space operator S n by Shimorin [3T] 
where it is called the Cauchy dual of the left invertible operator S n . For the latest 
developments, see [14] and the references there. 

Just as in the case of unshifted observability operators, it is possible to represent 
the range of a fc-shifted observability operator as a reproducing kernel Hilbert space. 

Theorem 5.6. The reproducing kernel for the space S k Ran (with inner 
product induced by A n (y) ) is given by 

Ak{z,Q = z k l h CR nih {zA)&-^ CtA R nik {C,AYC*. (5.17) 

Proof. Observe that by P~2T|) . 

OO 

S*D n , k ,c,AX = (T 1 ) {CAi- k x)zK 
J= k 

It then follows from (jTTTTjl and (|4~28l) that 

\\S k n D n , k ,c,Ax\\ 2 Aniy) = ( £ (A*y- k C*CA^ k x, x\ 

\o=k I x 

n+ j'+^ 1 ) (A*yc*CA J x, x\ 

3=0 I x 

= (®n,k,C,AX, x) X ■ (5.18) 

Therefore, by the general principle from [9] (see the proof of part (3a) of Theorem 
15.11 above), it follows that the reproducing kernel for S k Ran D n ,k,c, A is given by 
(|5.17|) as asserted. □ 

With all this model machinery in hand, we are now able to give the additional 
properties concerning rt-observability and exact n-observability. 

Proposition 5.7. Let (C, A) be an n-output stable pair. 

(1) // (C, A) is n-observable, then (C, A) is also k-observable for 1 < ft < n. 

(2) It can happen that (C, A) is exactly n-observable but there is a k with 1 < 
k < n such that (C, A) is not exactly k-observable. 

(3) // (C, A) is exactly n-observable (so <£> n ,o.c.A = Gn.c.A is strictly positive 
definite), then <S5 n ,k,c,A is strictly positive definite for k = 1, 2, 3, . . . . 

Proof of (1): It is not difficult to check by making use of the formula (15. ip that the 
adjoint Bergman shift S* has no nontrivial isometric vectors for n > 2: if n > 2 
and f G A n (y) satisfies \\S^f\\ — \\f\\, then f — 0. Making use of the intertwining 
relation (|5.5j) and the weighted Stein identity (|4.10j) . one can show: for n > 2, 
if (C, A) is n-observable, then (C, A) is also (n — 1)- observable. Then a simple 
induction argument leads to statement (1) in the proposition. 



' oo 

EC 
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Proof of (2): We let E: A n {y) ->■ y be the evaluation-at-zero map E: f i-> /(0). 
The pair (£7,5*) is n-observable since the kernel of O n ,E,s* is trivial and it is 
exactly n-observable since Gn,E,s* — lA„(y) 1S strictly positive definite. However, 
we have from (|4. 10[) 

Qn-1,E,S* = Gn,E,S* — S n Gn,E,S*Sn = I — S n S^ 

so that, according to (|5.1j) . G n -i,E,s* '■ z 3 \-> nXj-i z 3 . Thus the (n — l)-gramian 
Gn-x,E,S" is n °t strictly positive definite so that the pair (E,S^) is not (n — In- 
exactly observable. 

Proof of (3): If (C, A) is exactly n-observable, then Corollary 15.21 tells us that 
(C,A) is similar to a model pair (C, A) := (i?^, 5*|x) for an 5*-invariant sub- 
space C A n (y). If (C, A) and (C,A) are similar output pairs, then it is easily 
verified that the associated gramians &„ t k,c,A = Rn,k{BA)[C*C] and <& n , g ^ = 

-Rn,fc(-S^)[C*C] are congruent. Thus to show that <£> n .k.c.A is strictly positive def- 
inite, it suffices to show that & t ^is strictly positive definite. For the model 

pair (C, A) one can check that 

<& n.k.,C.A = P M^n,k,E,St\M- 

Thus strict positive-definiteness of <3 n k q ^ follows from strict positive-definiteness 
of & n ,k,E,s»- Observe that in view of (|4.27[) and (|5.4p . we have 

oo oo 

Therefore, 

oo 



3=0 3=0 r 



cfe»-» f \ /fwij -f ,?'+fc 

j=0 ^ n ^+ k 



and by computation (|5.18l) . we have 



from which we conclude that 



,'=0 ^ n '3+k 



®n,k,E,S* : £ /i^ "> £ Jt ^Si^- (5-20) 

As u„ 7 - is decreasing in j, we see that M "' J > 1 and hence ©„ k E s* ~> Ia (V) is 
strictly positive definite. □ 

6. Functions 0„.fc and metric constraints 

In this section we take a closer look at the transfer functions 0„.fc introduced 
by the realization formula (|3.8p . As was mentioned in Section 3, for an n-output 
stable pair (C, A), the associated backward-shifted observability operators O n ,k,c,A 
are bounded for all k > as operators from X into A n (y). In this case, the 
multiplication operator Me„ t given (according to f|3 . 8j) ^1 by 

Afe*.* - ( fc+ fc _1 ) A + S n D n , k+1 , c ,AB k : U k -> AiQO (6-1) 
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is also bounded. Therefore, the output function y in (|3.6|) . 

y(z) = O nt c,AX + y ] z k O n , k (z)u k 



N 

£• 

fc=0 



(6.2) 



belongs to A n (y) for every choice of x £ X and u k £ U k for each N = 1,2, ... . We 
next impose some additional metric relations on [ ^ ] > specifically one or more 
of the relations 



A*& nik+1}C}A B k + ( ) • C*D k = 0, (6.3) 

B* k <8 n , k+ i, c ,AB k + ( "+ fe fe - 1 ) • D* k D k < I Uk , (6.4) 

B* k <8 n , k+1 , c , A B k + (™ + £- x ) ■ D* k D k = I Uk , (6.5) 

and show how these lead to boundedness and orthogonality properties for the as- 
sociated multiplication operator Me fc . Due to equality (|4.30|) . it turns out that 
relations (|6.3|) and (|6.4[) are equivalent to the matrix inequality 



'A* 


C*~ 




B* k 







<&n,k+l,C,A 









( n r)'^. 



A B k 
C D k 



< 



<$n,k,C,A 









while the equalities 



and (16.51) are equivalent to the matrix equality 



'A* 


C*~ 




B* k 


D l_ 





®n,k+l,C,A 









A B k 
C D k 



<$n,k,C,A 









(6.6) 



(6.7) 



The two latter conditions are of metric nature; they express the contractivity or 
isometric property of the colligation operator \J k — [ ^ p k ] with respect to certain 
weights. Note that in the classical case n = 1, the formula (|6.10l) below amounts 
to the following well-known fact: if the colligation operator U 
and if we set Q(z) = D + zC(I - zA)~ 1 B, then 

z-e(c)*e(z) 



— [cd~\ i s isometric 



1-zC 



= B*{I-CA*y 1 {I-zA)- 1 B. 



(6.8) 



Lemma 6.1. Let (C, A) be an n-output stable pair and let <d n _ k be defined as in 
for some integer k > and operators B k £ C(lA k , X) and D k £ C(U k , y). 

(1) If equality (|6.3I) holds, then 

(a) O n> c,A x is orthogonal to S^Q n ^ k u for all x £ X and u £lA k . 

(b) S k Q nik u is orthogonal to S™Q n ^ k u' for all m > k and u,v! £ lA k . 

(2) Moreover: 

(a) If inequality (|6.4[) holds, then the operator S^MQ n k is a contraction from 
U k into A n (y). 

(b) If both (|6.3[) and (|6.4I) hold, i.e., if (|6.6|) holds, then the operator S^MQ n k 
is a contraction from H 2 (U k ) into A n (y). 

(3) Similarly: 

(a) If equality (I6.5P holds, then the operator S*M@ n k is an isometry from U k 
into An{y). 

(b) If (JO]) and dfO| hold, i.e., if ([6/7]) holds, then 



\Sn®n,kf\\\ n (y) 



II/" 2 



HHu k) - E IK 7 " S* n S n ) l ' 2 S k n e n>k S^ff (6.9) 



for every f £ H 2 (U k ). 
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(c) If {EH) holds, then 

- e n , fe (z)*e n ,fe(C) =^ k BiR nM {zAy<8 nM+lfiA R nM (CA)B k (6.10) 

f^>n,k 

- z^^ n ^BlR n . k+ i(zA)*<8 n ^ kfi ,AR n ,k+i{C,A)B k . 
Proof of (1): We first observe the power series expansion 

oo 

©n,feW - ^\D k + M;, - +fc+1 • CA^*^ 1 . (6.11) 

i=o 

which is an immediate consequence of formulas (|3.8j) and (|4.2T[) . We then make use 
of expansions (|3.7[) , (|6. 1 1|) and the definition of the inner product in A n (y) to get 

(S^Q nik U, O n ,c,Ax) An{y) = (J>n,k ■ {ll~^ k D k U, fX^ k CA k xj^ 

OO 

+ J2»n,j+k+i ■ U-) +k+1 CAiB k u, H-} J+k+1 CAi+ k+1 x\ 
j=o 

•n+*-i) C *D k +A* (j2 {l+fXi) A* j C*CA^j B)j u, A k ^ 

= ({{ n+k k - 1 )C*D k + A*& ntk+lt c, A Bk)u,A k x) x (by ^M) 
= (by flO])). 

which proves part (a). Verification of part (b) is quite similar: for m > k we have 
(SZ l e n , k u', S k Q n , k u) An(y) = ti n , m ■ (tf k D k u\ fi^ m CA m - k - 1 B k u) y 

+ £/*»j+m+i • {^ +k+1 CA^B k u', ^ j+m+1 CA*+™- k - l B k u) 

(n+fc-i) c*D k + A* \ ^p(jXlX\) A* j C*CA^j B k J u', A^^B^ 
(((n+fc-i ) C *.D fe + A*C n)fc+1)C) A5 fc ) u', A m - k - l B k u) x = 0. 

□ 



Proof of (2): According to (jOj) , 

( fc+ r 1 ) • Pfc«lly + (Cn.k+wM < Nlw» (6-12) 

for all u &U k . We now may make use of (|5.18|) to get 

\\Sn®n,kU\?A n{ y) = IK( fe+ r X ) + S k+1 D n , k+1 , C ,ABk)u\\ 2 An{y) 

= II ( fe+ r X ) S k n D k u\\ An{y) + \\S k+1 Q n . k+1 .c. A B k u\\ An{y) 

= \\D k u\\y + {® n , k+1 , c ,A B kU, B k u) x < \\u\\l h . (6.13) 

Thus, S k M@ n k is a contraction from U k to Ai(J^) and part (a) follows. To prove 
part (b) , following ideas from [29] , under the assumption that both (|6.3j) and (|6.4j) 
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hold, we shall show that for any Z/4-valued polynomial f(z) = fjz\ 

3=0 

m 

\\s k n Qn, k f\\ 2 Aniy) < \\f\\ 2 m ( u k) = E WfiWk- ( 6 - 14 ) 

3=0 

Let S* be the operator of backward shift on H 2 (Uk) so that for the polynomial / 
m— 1 

as above, (Slf)(z) = fj+iz^. By statements (lb) and (3a) of the lemma, we 

j=o 

have 

2 



|<S^6n,fc/|| — 



3=0 



<ll/o| 



3=1 
m — 1 



(by (lb)) 



3=0 

s* +1 en, k stf\\ 2 

2 



-ll/oll 2 

= ||/o|| 2 + ||^e n , fe ,s 1 */| 

Replacing / by S\ J f in (|6.15[) gives 

||5X,^rvil 2 <ll/,-|| 2 + ll^e n , fc (s 1 *y+ 1 /ll 2 

-H(/-s:s„)i^e n ,fcW +1 /ll 2 

for j = 1, . . . , m. Iteration of the inequality f|6.15D using (|6.16|1 then gives 

m m 
3=0 j=l 



(6.15) 

(6.16) 
(6.17) 



< 



3=0 



Letting m — > oo in (|6.17[) now implies the validity of (16. 14|) for every / £ H 2 (Uk) 
and the proof of part (b) of (2) is now complete. □ 

Proof of (3): In case (I6.5[) holds, then (|6 . 1 2[) holds with equality and part (a) of 
(3) follows. If also (|6.10[) holds, then f|6 . 1 5[) holds with equality as well in (|6.14l) . 
(j63B]) . and (|6TT7| . Part (b) of (3) now follows by letting m oo in (|6TT7) . 

It remains to verify the formula f|6. 10[) under assumption (|6.7I) . The identity 
(16. 7[) is equivalent to the collection of identities 

0*0 C*Dk 

A*<8 n ,k+1,C,AA -\ = ®n,k,C,A, A* & n ^+l,C,AB k H = 0, 



Bl<$ n ,k+l,C,ABk 



Hn,k 



L U k - 



(6.18) 
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We use these relations to compute 



f^n,k 



e„, fc (z)*e„, fe (c) = 

+ zB* k R 7l!k+1 (zA)*C* 

f^n^k 

D* k D k 



Dk 



(CR ntk+1 ((A)B k 



2 zBlR n , k+1 (zA) 

f^n,k t^n k f^n.k 

-zCB* k R n , k+1 (zA)*C*CR n , k+1 {(A)B k 
1 



*°* Dk t^R n .k + i{QA)B k 



t^n,k 



-B k & n ,k+l,C,ABk + zBlR n . k+ i(zA)* & n ,k+l,C,ABk 
l^n,k 

+ (Bl<3 n . k+ i y c,AAR n ,k+i((A)B k 

- (J, nt kZCB k Rn,k+l(zA)* (<& n ,k,c,A ~ A* ® n ,k+l,C,A A ) R n ,k+l{C A ) B k 

where we made use of (|6 . 18[) in the last step. By making use next of relation (j2 . 14[) 
we can continue the computation as 



f-Ln,k 



e n , fc (.z)*e n , fc (c) 



-B* k & n . k+ltC ,ABk + B* k R n . k {zA) 



f-i>n,k 



n,fc+l,C,A-Bfc 



+ B* k <3 n , k+hC ,A RnA< A ) ~ 



+ Vn,kB k R nk (zA) 



I 



fJ-n,k 



Bl 



^n,k 



®n,k+l,C,A Rn,k((A) 



B, 



~ Hn,kZ(B k R n , k+1 (zA)*G n , k ,c,ARn,k+l((A)B k 
— Hn,kB k R n , k (zA)*<3 rlik+ i i C,ARn,k{(A)B k 

— ^n,kzC,B k R n ^ k+ i(zA)*Q n ^ k ^c,ARn 1 k+l{C A )Bk 

verifying formula (|6.10j) . 

The following result is an immediate consequence of Lemma 16.11 



□ 



Corollary 6.2. Let us assume that the pair (C, A) is n-output stable and that 
relations (|6 . 3[) . (|6.4|) hold for all k > 0. Then the representation (|6.2p of the 

function y is orthogonal in the metric of A„ (y) and 

oo oo 

\\y\\A n (y) = \\0„,c,ax\\ 2 + J2\\ e ^ u *f ^ IKc.axWx + T, H u *H"* • ^ 6 - 19 ) 



fc=0 



// relations (|6.4[) hold with equalities for all k > 0, then equality holds in (|6.19[) . 

Observe that in case the pair (C, A) is exactly n-observable (so that <3 n ,k,c,A is 
strictly positive definite for all k > by Proposition 15. 7|) . the inequality (|6.6j) can 
equivalently be expressed as ||S|| < 1 where S is the operator given by 



<S 1/2 






'A 


B k ' 






C 







®nlk, 2 C,A 

Iu h 





X 




X 




u k _ 


— > 


y 



(6.20) 
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Another equivalent condition is that ||^*| < 1 which in turn can be expressed as 



'A 






C 







n,k,C,A 








In 





'A* 


C*~ 


< 









Hn,kly, 



3.21) 



Note that equality (|6.7j) means that the operator S is isometric. Of particular 
interest is the case where S is coisometric, i.e., where the colligation operator U k = 
[ C D k ] i s coisometric with respect to the weights indicated below: 

-l 



'A 


B k ~ 




C 







n,k,C,A 

Iu k 





'A* 


C*~ 














Hn,k,Iy_ 



3.22) 



We can now derive the shifted weighted Bergman-space analogue of the familiar 
identity for the de Branges-Rovnyak kernel which is a dual version of (|6 ,8|) : if 
U = [q £,] is coisometric and if we set Q(z) = D + zC(I — zA)~ 1 B, then 

J-9(z)e(0* = c ,j _ zA) -i,j _ £ A *yi c *, 
1 - z( 

Lemma 6.3. Let (C, A) be an exactly n-observable n-output stable pair and let 
Q n ,k be defined as in (|3.8p for some operators B k E C(l4k,X) and D k E C(U k ) 
subject to equality (|6.22|) . Then 



Iy 

l^n,k 



e» )fc (z)e n , fe (C)* =GR n , k {zA)&-l^ A R n>k ^A)*C* 



- z( ■ CR n , k+1 (zA)&-^ +hCjA R n , k+1 (CA)*C*. (6.23) 

Proof. The proof parallels the verification of the identity (I6.10|) done above. The 
weighted-coisometry condition (|6.22[) gives us the set of equations 

A®-lc,A A * + B * B * = C n,fc+l,C,X. C®n%,C,A A * + D * B * = 

C ®n,k,C,A C * + D kD* k = fj, n ,kly- 
We then compute: 

h 



3.24) 



f-t"n,k 



e„, fc (z)e„, fc (cr 



Iy 








Iy DkD* k 


f-t>n,k 





+ zCR n ,k+i {zA)B k 
- zCR n , k+1 {zA) 



D* 
u k 



(B* k R n , k+1 ((A)*C* 



-c- 



f-Ln,k f^n 7 k 

zC ■ CR n , k+1 {zA)B k BtR^ k+1 (C,A)*C* (by §M) 
-l 



c 



c* 



J n,k,C,A ,, 
f*n,k f^n,k 



C 

zCR^ k+1 {zA)A&~l kCA 



fJ>n,k 



+ <—®nl.C.AA*Rn,k + l(CA)*C* 



= c 



'n,k,C,A x 

t^n,k 

z( ■ CRn tk+ i(zA) 
Ix 



R n , k+1 {CA)*C* 



[J>n,k 



+ zR n , k+1 (zA)A <S n 



,,k,C,A 



Ix 

[^n,k 



(A*R ntk+1 ((Ay c 



- z( ■ CR n , k+1 (zA)&-' k+1C A R n , k+1 ((;AyC* 
CR n MzA)®-^ c%A R n M.(AyC* 
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z(-CR nik+1 (zA)<5 



-l 

n,k+l,C,A 



Rn, k+1 (CArC* (by (EH). 



□ 



Remark 6.4. More generally, if O n ,k(z) is given by (|3.8[) and if we do not assume 
the weighted coisometry condition (|6.22j) . then the decomposition (|6.23j) holds in 
the more general form 

Iy 



H-n,k 



@n,k{z)®n,k(0* = CRn,k(.zA)<S;Xc,A R ^A)*C* 



-z(-CR n ,k + i(zA)<8-i k+1 
where the defect kernel £) is given by 



CA R n , k+ i((Ayc* + ~ k (z,C) 



^k(z,() 



zCR n ,k(zA) 



iy 



<3 



n,fc+l 








'A 


Bk 




'«»,* o" 




'A* 


C*~ 




C 


D k 




Iy 




Bt 




)■ 



(^((Ayc* 

Iy 

Hn,k 



Since equality (|6.22[) implies inequality (|6.6p . it follows that under assumption 
of Lemma T6.31 all the conclusions of parts (1) and (2) in Lemma [6.11 are true. To 
have all conclusions true, we need the operator (|6.20l) to be unitary. 

Lemma 6.5. Suppose that we are given an integer k > and an exactly n- 
observable n- output- stable pair (C, A) with A £ C(X) and C S C(X,y) . Then 
there exist operators B k G C(Uk, X) and D k 6 C(Uk,y) such that equalities (|6.22[) 
and (|6.6p hold. Explicitly, such B k and Ck are essentially uniquely determined by 
solving the Cholesky factorization problem: 



Bk 



[Bl D* k ] 



<3 



n,k+l,C,A 







A 




C 



(55 



n,k,C,A 



[A* C*] (6.25) 





Hn,kly_ 

subject to the additional constraint that the coefficient space Uk be chosen so that 
[ D k ] : Uk — > X © y is infective. 



Proof. By Proposition @TTTJ the weighted Stein identity (I4.30[) holds for each k > 1. 
Since (C, A) is exactly observable, the gramian <S n .k.c.A is strictly positive definite 
and then it follows from (I4.30|) that the operator 



r l 

<5 



_ l ' _i 



x -> x®y 



is an isometry. By extending this operator to a coisometric operator (I6.20[) we 
arrive at B k and D k solving (|6 . 25[) . Further, extension of this operator to a unitary 
amounts to the additional restriction that [ D k k ] be injective. □ 

7. Beurling-Lax theorem for A n 

We now present our four approaches to Beurling-Lax representations for the 
Bergman space setting; these were sketched for the classical case in the Introduction 
above. 

We start with a general observation. If the subspace M C A n {y) is ^-invariant, 
then is ^-invariant; hence we may apply Theorem 15.11 part (3) (with in 
place of M viewed as sitting isometrically inside A n [y )) to conclude that there is an 
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rt-isometric pair (C, A) so that Ai^ — R&nO n ,c,A', in fact, as indicated in Theorem 
5.11 we may take (C, A) to be the model output pair (C, A) = (E\ M ± , S^\ M ±), and 
O n .c.A amounts to the inclusion map of Ai x into A n (y)- Since ||0 n ,c,A2:||^ r y \ = 
(Gn,C,A%, %)x for every x G X, it follows that Ai 1 ' is a reproducing kernel Hilbert 
space with reproducing kernel 

k M ± (z, o = c(i zA)- n g^ c ji - U*r n c*. 



It then follows that Ai = (Ai ± ) ± has reproducing kernel 



kjvi{ z X) 



(l-zO n 



- C(I - zA)- n g-'(i - (A*)~ n C 



(7.1) 



7.1. The first approach: partially isometric multipliers. Suppose that we 
are given a shift- invariant subspace M contained in the Hardy space H 2 (y). As 
we just observed in the previous paragraph, then Ai 1 - is backward-shift-invariant 
and hence can be represented as the range Ran Oc,a of an observable output pair 
(C, A); we can even take (C, A) to be the model output pair: C: f h-> /(0) and A 
equal to the backward shift restricted to X = M. , From the identity 

Gi,c,A-A*g hC>A A = C*C (7.2) 
(where Oi.c.A is the identity operator in the model case), we see that the column 

' ! AC-V 2 1 

is isometric. Hence there exist operators Bi : U\ — > X 



matrix 



G 1/2 AG~ 1/2 



(= AA ) and Di : U\ — » y (where U\ is an appropriate coefhcient Hilbert space) so 
that the operator matrix 



1 




X 




X 






u 


— > 


y 



(7.3) 



@i,c,aAGi C,A 

_ 1 
Cy\,c,A 

is unitary. We then set 

*i(z) =Di + zC(I - zA^Bl (7.4) 

Making use of the unitary property of the matrix (|7.3[) . we deduce that is 
inner in the sense that its boundary values exist almost everywhere on T and are 
coisometric, and that the kernel identity 

Iy 



1-zC 



= c<(i-zA)- 1 g- 1 CA (i-(A*)- 1 c* 



(7.5) 



holds. For the case where n = 1, the latter identity implies the factorization of the 
kernel (fTTT]) 

*i(*)*i(0* 



kM(z,0 



(7.6) 



1-zC 

Standard reproducing-kernel arguments (details appear in the proof of Theorem l7.ll 
below for the more general representation (I7.8P ) then imply that Ai = -Ai(U) = 
• H 2 (U), i.e., serves as a Beurling-Lax representer for Ai. One can view this 
approach as a state-space adaptation of the de Branges-Rovnyak reproducing kernel 
approach (see |12[ 113)): it has proved to work well in a number of multivariable 
settings (see e.g. 018]). The next theorem presents a Bergman-space analogue of 
the representation (|7.6[> . 
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Theorem 7.1. Let M. be a closed S n -invariant subspace of A n (y)- Then there is 
a partially isometric multiplier 



M F =[M Fl ■■■ M Fn ] : Aj {Uj ) A„ (y) 

3=1 

(so F(z) has the form F(z) = \F\(z) ■ • ■ F n {z)\ with each Fj holomorphic on B 
with values in jC(Vij, y) for appropriate coefficient Hilbert spaces Uj, (j = 1, ... ,n), 
so that 



M = M F l®A j (U j ) 



(7.7) 



J =1 



Proof. To prove the theorem, it suffices to produce coefficient Hilbert spaces U\ , . . . , U n 
and holomorphic operator- valued functions F 3 w D — >• £(Uj, y) so that 



(7.8) 



where kj^ is given in (|7.1[) . Indeed, assume that (|7.8|) holds. Recall that Aj is the 
reproducing kernel Hilbert space with reproducing kernel kj(z, Q = (1 — zC) - ' 7 an( i 
define the operator X : M = Ti{k,M) — > ©™ = i AiQii) by 



'ki 0^1(0*1/" 

k„(-,C)^n(C)*y 



(y e y, c g 



Then we have for yi^y and £j € D, 



/ iV 

E 



k,(-:0)/ r i(0)V 



k„ (-,0)^(0*% 



>E 



'ki(-,C0^i(Ci)*W 

.kn(-,Ci)^n(Ci)*W 



AT / AT 



E (E^COMCi, 0)^(0)*%, w 



A 7 " 



E ( k M(0,(j)yj, v*)y = 

i,3=l 



N 



E^("'0)% 

3=1 



We conclude that X extends uniquely from the span of the kernel functions fc» (•, £)y 
(C 6 D, y e ^) to dchne an isometry from %{kM) = M into the direct sum space 
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0"=i Ai{Ui). Moreover, we compute, for hi £ Ai(Ui) (£ = 1, . . . 



X* 



(0,2/ =(X* 



hi 



kM(;,Qv 



xk M (-,()y 



ki(-,C)^i(C)*y 



)" =1 ^«(^) 



from which we conclude that 



hl(z) 



= (F 1 (C)h 1 (0 + --- + F n (0h n (0,y) y 
^ Fi{z)hi(z) + ■■■ + F n (z)h„(z) 



X* = Mp : 

U„(*)J 

is the multiplication operator Mp = \M p 1 ■ ■ ■ Mp n ] induced by the matrix 
function F(z) — \Fi(z) ■ ■ ■ F n {z)\ and Theorem 1 7 . 1 1 then follows. 

Thus the proof of Theorem 17.11 is reduced to the construction of holomorphic 
operator-valued functions Fi(z), . . . , F n (z) so that (|7.8p holds. The construction 
proceeds via an iterated unitary completion process as follows. 

For j — 2, . . . , n we do a construction similar to that in (|7.2[) - (|7.4p but based on 
the identity (see P~TD|) ) 

Qj,c,A - A*Q^ c ,aA = Gj-\,c,A (j = 2, . . . , n). 

We find operators Bj : U n +\-j — > A" and Dj : U n +i-j — > y so that 

" ^ B/ 



^C,A 



and 



B* D* 



8j,C,A 
Jy 



B* 1 
1 B, 

1 





^•,C,A 


" 









Jj,C,A 












(7.9) 
(7.10) 



(7.11) 



In fact, the latter equalities determine Bj and Dj uniquely up to a common unitary 
factor on the right: 

■5 = iQjk A M;J-. a A-}>. 1) = -Q;j 1 ,c,A^Gj,c,AB j . (7.12) 
We now define the functions 

*,(z) = D, + zg]_ l c A {I - zA)- l Bj (7.13) 
for j = 2, . . . , n, which are inner and satisfy the identities 

— -= = Qj-^cAK 1 ~ zA > j,c,A\ - Q A ) y/-l,C,A- (<- 14 ) 

We finally define operator- valued functions Fg(z) : Ui ~^y via 

F < («) = C(J-2!A)-( n -')g~i )C)A * n+ i_ < («) for ^ = l,2,...,n-l, 
F n (z) = *!(z) = D x + zC(I - zA)- x B! (7.15) 
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and claim that this choice of F± , . . . , F n satisfies (|7.8[) (with coefficient spaces Uj — 
X for j = 1, . . . , n — 1). Indeed, substituting equality (|7.14|) (for j = n) into the 
formula (|7.ip for fc» and then making use of formula (|7.15l) for F%, we have 



kM(z,C) 



h 



(!-<)' 



(l-zC) n 



c(i - zA)-^g n \ 



C.A 



j-* w (*)*»(0* 



1-zC 



1-zC 



Q-l 

C(I - zA)- {n - x) "-^ (j - (A*)- {n -^C* 
1 - z( 



One can now use an inductive argument to show that in general 

k 



/ >n ly FA 

(z, C) = : ^-=— + > — 

1 S; (l-zC) n f-f (1 



WW 



t (l-*cy 



(1 - zC) 



(7.16) 



for k = 1, 2, . . . , n — 1. For the final step, we note that the last term in (|7.16[) with 
k = n — 1 is given by 



Q-l 1 

C(I-zA)- 1 ^ — (J-CA*)- 1 ^ = = 

(1-zC)"" 1 (1-zC)™- 1 



7y - * 1 (z)* 1 (C)' 



l-*c 

Iy_ F n {z)F n (Q* 

(1 - zC)™ (1 ~ <)" 

Plugging this expression into the fc = n — 1 case of (|7.16l) leaves us with 
wanted. 



as 
□ 



There is an alternate formula for vf^- (and hence also for Fi) which may prove 
useful in applications. 



Proposition 7.2. The function tyj in (|7. 1 3[) for j = 2, . . . , n can alternatively by 
given by 

= - zA)~%\zI A*)BJ\ 

Hence Fg(z) can alternatively be given for t = 1, 2, . . . , n — 1 fry 



(7.17) 



F,(*) = C(I - zA)-^gj e c A y n+1 ^(z) 

= c(i - zAy^-^g-l^izi - A*)B^ X _,. 



(7.18) 



Proof. In the proof, we shorten notation Qc,A,j to (/,*. We substitute the second 
equality from (|7TT2"|) into (|7TT5|) to get 

= - gjj^GjBj + zg]_ x {i - zA)-^ 

=gj_ 1 (i - zA)- 1 [-(/ - zA)gj\A*g 3 + */] b, 

=$*_!(/ - zA)- 1 [ z (i + Agr} lA *g 3 ) - g-\A*g,\ b,, (7.19) 
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We next observe the equalities 

Qj\A*Qj (I + AGr^Gj)- 1 = Gi 1 A*g 1 , 

and 

g j + QiAQj^A'Qi = (g^ - Ag^A*)- 1 , 

which both follow from (|7.9|) . We then use these equalities along with (|7.19[) and 
the first formula in ([7715} to get ([77171) : 

=Q]-i{i - zA)- 1 [zi - gj 1 A*g j \ (i + Agr^g^Bj 
=g}_ 1 {i - Z A)-%\zi - a*){q s + c;,.u; ; \.\g,:B, 
=gj_ 1 (i zA)- x gT\zi - A^gj 1 - Agj 1 A*)~ 1 B j 
=gf_ 1 (i - zA)- l gj\zi - A*)Bj\ 

□ 

Given a shift invariant subspace Ai, let us refer to the function F — \Fi ■ ■ ■ i 7 ^] 
constructed as in Theorem 17.11 as the partially isometric Bergman inner function 
associated with Ai. 

7.2. The second approach: shift-invariant subspaces contractively in- 
cluded in A n (y). Let us say that the Hilbert space Ai is contractively included 
in the Hilbert space H if Ai C H as sets and moreover ||m||-^ < ||to||wi for all 
m G Ai. We also say that an £(U , y)-vahied function is a contractive multiplier 
if the operator Me : f(z) h-> 9(z) • f(z) of multiplication by O defines a contractive 
operator from A n {U) to A n {y). 

Theorem 7.3. A Hilbert space Ai is such that 

(1) Ai is contractively included in A n {y), i.e. Ai C A n (y) with ||/||^ ^ < 
\\f\\M for all feM, 

(2) Ai is S n -invariant: S n Ai C Ai, 

(3) the operator A = (S n \M)* is a strongly stable n-hypercontraction, 

if and only if there is a Hilbert space IA and a contractive multiplier so that 

M = e-A n (U) (7.20) 

with lifted norm 

lie • f\\ M = WQfW^u) (7.2i) 

where Q is the orthogonal projection onto (Ker Mq) 1 - . 

Proof. We first verify sufficiency. Suppose that Ai has the form (|7.20p for a con- 
tractive multiplier O with .M-norm given by (|7.21[) . Since ||Me|| < 1 it follows 
that 

ll©/IU„w = WeQf\\ An (y) < IIQ/IU„(w) = l|e/||M 

i.e., (1) holds. From the intertwining equality S u Mq — M<s>S n , property (2) follows. 
The latter intertwining equality also implies MeSVilKcr A/ e — which can be written 
equivalently in terms of the orthogonal projection Q onto (Ker Me) 1 C A n {U) as 
QS n (I - Q) = 0. Thus, we have 

QS n = QS n Q and 5*Q = QS*Q. (7.22) 
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Furthermore, for every f,gE A n (pi), we have 

(Qg, AQf) M =(S n Qg, Qf) M = (QS n g, Qf) M = (QS„g, f) An (u) 

=(QS n Qg, f) An (u) = (Qg, SlQ^Aniu) = (Qg, @S*Qf) M , 

which implies that A : Of — » QS^Qf. Iterating the latter formula gives 

A j : 6/ -> eS*JQf for j > 0. (7.23) 

Since S* is a strongly stable n-hypercontraction on A n (y), we conclude from (|7.22p . 
(|7T23)) . and (|7T2T|) that 

lim ||^6/|| M = lim \\QS* j Qf\\ M 

iV— >-oo iv— »oc 



: lim \\QS?Qf\U(u) = I™ ||^Q/|U n(W ) = 



and also 



for A; = 1 , . . . , n which shows that A is strongly stable n-hypercontraction on M. 
and therefore completes the proof of sufficiency. 

Suppose now that the Hilbert space M. satisfies conditions (1), (2), (3) in the 
statement of the theorem. Using hypothesis (2) we can define the operator A := 
(S n \M)* on M. and since it is n-hypercontractive by hypothesis (3), the operator 
YA,n(I) is positive semidefinite. Choose the coefficient Hilbert space U so that 

dimbf = rankr j 4.„(J) 

and then choose the operator C : M — > U so that 

C*C = T n>A (I) = Y,(-iy (f) A* j A* > 0. 

3=0 W 

Then (C, A) is an n- isometric pair and, since A is strongly stable by hypothesis (3). 
it follows from part (2) of Proposition 14. 71 that the observability operator 

O n , c ,A-- f^C(I-zA)- n f 

is an isometry from Ad into A n (y). By part (1) of Theorem I5.1[ we have the 
intertwining equality S^O rh c,A = O n ,c,AA. Taking adjoints then gives 

(O n ,c,A)*S n = A*(O n ,c,A)*. (7.24) 

The inclusion map i: M, — > A n (y) is a contraction by hypothesis (1). Moreover, 
i A* = S n t : M -> A n (y). Therefore the operator 

r = i o (o c ,A,nT ■ Mu) -> My) 

is a contraction and 

RS n = i(O n , c ,A)*S n = iA*(O n , c ,A)* = S n L(O n , c ,A)* = S„R. 

Therefore (see [28 ) R is a multiplication operator, i.e., there is a contractive multi- 
plier so that R — M&. Since O n .c.A ■ M — > A n (y) is an isometry, it follows that 
Ran(O n! c,A)* = M and also that M = 6 • A n {U) with .M-norm given by (|7.21|) . 
This completes the proof. □ 
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Remark 7.4. It is of interest to consider Theorem l7.3l for the case where condition 
(1) is strengthened to 

(1') A4 is isometrically contained in A n {y). 

For the special case where n = 1 (so Ai(y) is the Hardy space H 2 (y)), it is not dif- 
ficult to see that condition (3) is a consequence of (1') and (2). Then the representer 
is a coisometric multiplier and the proof of Theorem 17.31 reduces essentially to 
that of Rosenblum-Rovnyak (see [3D]) for the case n = 1. The same approach has 
been adapted to the multivariable setting of the Drury- Arveson space in Arveson [5] 
and McCullough- Trent [24] as well as the freely noncommutative lift of the Drury- 
Arveson spaces to the Fock space Theorem 2.14]. For the case n > 1, there 
is a new phenomenon: it is not the case that (1') and (2) imply (3): indeed, it is 
known that any Hilbert space operator can be recovered up to unitary equivalence 
as the compression of the Bergman shift S2 to the orthogonal difference Ai Q Af of 
two nested ^-invariant subspaces Af C Ai; in particular, the 2-hypercontractivity 
property of <S£ is n °t preserved when one considers (S^Ia^O* for an ^-invariant sub- 
space Ai. A second manifestation of the inapplicability of Theorem 17.31 to general 
isometrically-included S^- invariant subspaces of A n {y) is given by Theorem 17.11 
S^- invariant subspaces of A n (y) are modeled by a coisometric multiplier acting on 
a direct-sum Bergman space Aj(Uj) rather than just on a single Bergman 

space An(U). 

7.3. The third approach: Bergman-inner families. Given an S^-invariant 
subspace Ai of A n (y), from the fact that 



one can see that we always have the orthogonal-sum decomposition for A4: 



The space S%A4 has reproducing kernel (z,() —> z k ( kM(z,()> but since the op- 
erator S k : A n {y) — > A n {y) is only a contraction and not an isometry if n > 1, 
the metric induced by the above reproducing kernel is different from the metric 
of A n {y). To get the reproducing kernel for S*A4 consistent with the metric of 
An(y), we need to consider backward-shifted version £> n ,k,c,A of the observabil- 
ity operator 0n,c,A and the backward-shifted version & n .k.c.A of the observability 
gramian Q n ,c,A introduced in Section|4]above. It turns out that the space (S^M) 1 - 
has a nice characterization in terms of the operator £) n ,k,c,A- 

Proposition 7.5. The space (S^AA) 1 - is characterized as 



Pi s^m c P| s k n A{y) = {0} 



kgeo fc>0 



M = (S k M e S k+1 M) . 



(7.25) 



(S k M) 




(7.26) 



where we identify the first term with the subspace of polynomials of degree at most 
k - 1 in A n {y). 
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oo 

Proof. We wish to characterize all functions f[z) — fjZ-* which are orthogonal 

3=0 

fc-1 

to S k M. in An(y). We may write f(z) — p(z) + z k f(z) where p{z) = fjZ^ and 

3=0 

f G A n (y)- Clearly polynomials of degree at most k — 1 are orthogonal to S k A4, 
so it suffices to characterize which functions of the form z k f(z) are orthogonal to 
S k M. To this end, observe that S k f is orthogonal to S k M if and only if the 
function (S k )*S k f belongs to M ± = Ran O^ca- It follows from the formulas 
(HEU) and (fTT2]) that 

OC OO 

(s k n ys k n -. X/-' — 

3=0 3=0 ' J 

We thus conclude that S k f is orthogonal to S k Ai if and only if there exists a vector 
x G X such that 

oo oo , 

E 7jZ j = C(I - zA)- n x = J2 ( — • CAj 

Equating the corresponding Taylor coefficients gives 

fj = —-— ■ CA'x = ( n+j A_ k k ~ l ) • CA j x for all j > 0, 
and therefore, 

00 00 

f(z) = fa j = E (( ) • CA ' X ) Zj = ®n,k,C,AX, 

3=0 j=a 

by (|4.27p . Thus, / £ RaniDn^cA- As the analysis is necessary and sufficient, the 
result follows. □ 

With characterization (|7.26[) in hand, it is straightforward to derive the kernel 
function for the space S k A4 with respect to the metric inherited from A n (y). 

Proposition 7.6. Let M. be a closed shift-invariant subspace of A n {y) with re- 
producing kernel kj^ given by (|7.ip . Then the reproducing kernel functions for the 
closures of S k Ai and of S k Ai QS k+1 A4 in the metric of A n {y) are given by 

(n /l+k-2\ j \ 
E (1 -l^n-ili CRnMzA^-^RnMCATC*] , 

(7.27) 

k s*Mes^ M ( z >0 =z k t (("T 1 ) • Iy - CR n , k (zA)&-^ c>A R ntk (CA)*C* 

+z(CR n , k+1 (zA)®-i k+hCA R n , k+1 (CAyC*) . (7.28) 



Proof. By Theorem 15 . 61 the kernel & for the space S k K&nO n , k ,c,A (with inner 
product induced by A n (y)) is given by (|5.17l) . From the formula (|7.26|) for (SnM.) 1 - 
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we deduce that 

S k n M = f Siyj f) (S* Ran£)„, fe ,c,A) X 

= S^OO e (5* RanD„ AC> A). (7.29) 
Since the reproducing kernel for the subspace S 3 n y of A n {y) is z^C* f ' 

A > . \ 

we deduce that S^A n {y) = AiCV) 9 S J n y has reproducing kernel 



fca^ooM = 77^77 - E ( n T x ) ^ 



i=0 



OO / oo 

= E ("T 1 ) ^» = - fe c fe E ("li"" 1 ) ^ 

j = fc ' \j=0 

If we now make use of the Chu-Vandermonde identity (|2.6I) . we can continue the 
computation as follows: 

(oo n \ 
EE(tr 2 )-(T)^ 'U 
j=o e=i J 

= A fe E(^ 2 ) E("T0^ J K 

= ^£ ( ;_yr (7-30) 

From the representation for S'k.M in (|7.29l) . we see that 

^mO^O = k S k An{y) {z,Q -Sl k (z,C) (7.31) 

where &k is the reproducing kernel for Ra,nO n ,k,c,A- Combining this with (|7.30[) 
and the formula (|5.1T[) for J^. leads to the formula (|7.2T[) for k S k M . 

It remains to verify the formula (|7.28[) . Toward this end, note that in view of 
(177311) and COp , 

= k s*A n (y)( z X) - k s* +1 A n (y)( z '0 ~ &k(z, C) + &k+i(z, 
= k s*y( z , C) - &k(z, C) + C) 

= Z k f ( ) • Iy - & k (z, C) + 0- 

Plugging (15.17P into this last expression leads to f|T. 28[) as wanted. □ 

Lemma 7.7. Given an integer k > 1 and an exactly observable n-output stable 
pair (C,A), construct operators Bk € £(Uk,X) and Dk £ C(Uk,y) as in Lemma 



WEIGHTED BERGMAN SPACES 



41 



and let Q n ,k be the associated function given by (|3.8j) . Then the kernel (|7.28|l 
can be factored as 

k s*Mes^M( z > = z k f&nAz)QnA(y- (7-32) 

Proof. By Lemma f6.3[ identity (|6.23j) holds. Multiplying both parts of (|6.23p by 
z k C, k and combining the obtained equality with (|7.28|) we easily get (|7.32[) . □ 

Let us say that a family of operator-valued functions {Q Uy k ■ © £(Uk, y)}tLo 
is a Bergman-inner family if, for each k = 0, 1, 2, . . . , we have: 

(1) M Bn k : S\U k -t A n (y) is isometric, 

(2) M Qn k (S k U k ) is orthogonal (in A n (y)) to span{Me n >t S{Ui : t > k}, and 

(3) S k+1 M enk U k c 5pSH{S£Me Bi4 Z4: * > fc}. 
If {@ n ,k}k>o is a Bergman-inner family and we set 

oo 

M = e n , fc SfZ4 C A* (30, (7.33) 

fc=0 

it then follows that the multiplication operator 

oo 

M e = [M Qo M Bl M &2 • • •] : S k U k -)• A* (3-0 

fc=0 

oo oo 

maps the time-varying Hardy space H 2 ({U k } k >a) '■= S^Uk (where u = z fe Ufc £ 

fc=0 fc=0 
oo 

H 2 {{Uk\) is assigned the Hardy-space norm ||u|| 2 = || u fe|| 2 ) unitarily onto the 

fc=0 

£ n -invariant subspace M. C ^l„(3 ; ). Putting all the pieces together, we arrive at 
the following converse of all these observations which amounts to our third analogue 
of the Beurling-Lax theorem for the Bergman space setting. 

Theorem 7.8. Let M. be a closed S n -invariant subspace of A n (y). Then there is a 
Bergman inner family {Q n ,k}k>o so that Meff 2 ({W}t>o) (with Mq as in (|7.33|) ). 

Furthermore, given the S n -invariant subspace M. C A n (y), the associated rep- 
resenting Bergman-inner family {0n.fc}fc>o can ^ e constructed according to the fol- 
lowing algorithm: 

(1) Set X = M 1 - and define A E £{X) and C E C(X, y) by 

A = S*\ M ±, C/ = /(0) for feM X . 

(2) Construct [ ^* ] by solving the Cholesky factorization problem (|6.25[) in 
Lemma \6.5\ 



(3) Set e n , k (z) - ttt + zCR n , k+1 (zA)B, 



k- 



Alternatively, Bergman-inner families {& n ,k}k>o can be constructed from an ar- 
bitrary strongly stable n-hypercontraction as follows. Let A E C{X) be any strongly 
stable n-hypercontraction on a Hilbert space X and choose C E C(X, y) with dense 
range so that 

C*C={I-B A ) n [I x ] 

(One such choice is C = D n<A : X -> y where D nA = (ELoC -1 )* (*) A * k A k ) 1/2 
and y — T>n t A '■= RanDn.^.J Construct B k , D k as in Step 2 above and then Q n ,k 
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as in Step 3 above. Then {Q n ,k}k>o is a Bergman-inner family and any Bergman- 
inner family arises in this way. 

Example 7.9. By way of an example, suppose that M. — {0} C A n {y), so X := 
M 1 - = A n (y). Then 

A = S*, C = E:f€A n (y)^f(0)ey. 

In the proof of part (3) of Proposition (see ()5.20[) ). we saw that 

oo oo 

tf,.,.r.s; : J2 U# "> E ( 7 - 34 ) 

j=0 3=0 r ,jn 

and hence 

oo oo 

W5 ■■ E "> E (7-35) 
One can then work out that 

& n,k+l,E,S* _ S n^n„E,S^ S n = °) Mn,fc-fy ~ E ®n,k,E,S* E * = °- 

Hence, solving for [ ^ ] via the Cholesky factorization problem (|6.25|) (with (C, A) = 
(E,S*)) leads to ] = [§] for all fc, as is to be expected from the conclusion of 
Theorem EH with All = {0}. 

7.4. The fourth approach: wandering-subspace Bergman-inner functions. 

Suppose that M. C A n (y) is £ n -invariant. The subspace 

£ := M G SM 

has the property that £ C M and £ _L SViVW. If it is the case that £ generates M 
in the sense that M. = spanjS^ : fc = 0, 1, 2, . . . }, then one says that M. has the 
wandering subspace property (with wandering subspace equal to £). For the classical 
Hardy-space case n = 1, one has that £ = ©o • Wo for an appropriate coefficient 
Hilbert space Uq and the inner function Oo constructed in (|7.4p . Furthermore, M. 
admits an orthogonal decomposition 

oo 

M = S{£ (7.36) 

and this decomposition coincides with that in (|7.25|) . 

If n > 1, the subspace S 3 n £ is orthogonal to £ for all j > 1, but it is not orthogonal 
to >!?™£ if 1 < m ^ j. Thus, the best one can expect is the equality 

oo 

M = V Si£ (7.37) 

which indeed holds if n = 2,3; see [2], [31], [32] (see also [33] for a different ap- 
proach). Letting k = 1 in formula (|7.28[) and recalling that C5 n ,o,c,yi = Qn,C,A, 
O n ,o,C,A = O n ,c,A (see (|4.29l) ) and [i n ,o — 1, we conclude that the reproducing 
kernel for the subspace £ equals 

fc £ (*,C) =Iy- C(I - zA)- n Q-% A {I x -lA*)~ n C* 



zQC j - ^ I C- 1 ^ j - (A*)-' j C*. (7.38) 
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Following [27l [29] we say that the function O is a wandering-subspace Bergman- 
inner function whenever 

(1) Mq : Uo — >• A n {y) is isometric, and 

(2) 9 • % is orthogonal to S n Q • U for l>\. 

In case MqUq = £ is the wandering subspace for the SVi-invariant subspace M. , then 
we have the Beurling-Lax-type representation of A4 as the closure of O ■ C[z] ®Wo, 
where C[z] denotes the algebra of polynomials with coefficients in C. Construction 
of a Beurling inner function for the wandering subspace £ of M. amounts to focusing 
on the first element 0o in the Bergman- inner family associated with M. . Specifying 
Lemma I7T71 for the case k = then leads to the following. 

Theorem 7.10. Given a strongly stable n- output-pair (C, A), then there exist oper- 
ators B £ C(U, X) and D G C(U, y) which solve the Cholesky factorization problem 



[B* D-] 
Moreover, if O is defined by 



®n\.C.A 


" 




A 





h_ 




C 



0(z) = D + zC - zA)- j B (7.40) 
j=i 

where [£] solves (|6.25[) (with k = 1), then: 

(1) The factorization kg{z,C > ) = Q(z)Q(£)* holds and therefore, the multiplica- 
tion operator Mq maps U onto £ unitarily. 

(2) The subspace £ is orthogonal to S„@k ' U for every k > 1. 

(3) Me is a contractive multiplier from the Hardy space H 2 (Uq) into the Bergman 
space A n (y) ■ Hence, if M. has the wandering subspace property, we con- 
clude that M has the Beurling-Lax-type representation 

M = A n (y)-closure of M e H 2 (U). 

We note that the last statement in the theorem is a consequence of part (2b) 
of Lemma 16.11 We also note that realization formula (|7.40[) appears in [27l [29] ) in 
slightly different terms. 

7.5. The case of zero-based shift-invariant subspaces with one zero. In 

this concluding subsection we illustrate all four approaches by applying them to 
the simplest shift-invariant (zero-based) subspace 

M = {/ 6 An : f(a) = 0} (7.41) 

of A n where a is a fixed point in B. If n = 1, then 

z ~ — (y, 

M = b a (z)Ai, where b a (z) = — . 

1 — za 

If n > 1 and a = 0, the latter representation still holds. Let us assume that a^O. 
The one-dimensional space M 1 - is spanned by the function k„(z, a) = (1 — za)~ n 
and thus M. 1 - = Ran O n> c.A with A = a and any C ^ 0. We choose C to make the 
gramian G n ,c,A equal one. Thus, 

A = a, C=(l-H 2 )*. (7.42) 

Then a simple calculation reveals that 

Gj,c,A = (1 - \a\ 2 ) n ~ j for j = 0,...,n. (7.43) 
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An induction argument making use of the weighted Stein identity (|4.30j) combined 
with the initial condition & n ,0,C,A = &n,C,A = 1 then leads us to the formula for 
the shifted gramians: 

We may also relate the numbers (I7.44[) to the function R n ^k defined in (12.31) : 

<Sn,fc,c,A = (1 - \a\ 2 ) n R n M\a\ 2 ). (7.45) 
Formula (17.11) takes the form 

v ' (i-zC) n (i-zs)"(i-"C) n 

The first approach: Substituting (T712|) and (fTlBI into (|7TT2)| gives 
; = (1 - Icf) 1 ^, Dj = -a for j = l,...,n. 
We now get from (|T.15[) 

-Fn(z) = -a + = 6 Q (z), 

1 — za 

and we get from (|7.18[) 

ft * = n -J+ _< = M*)' H y~ fOT *=l,...,n-l. 

The identity (17.81) now amounts to 

i (i-MT 



&m(2,C) = 



{l-zQ n (1 - za)"(l - aC) r 
(1 - \a\ 2 Y b a {z)bJX) 



- (l-z5)'(l-a0 J (l-«0 n_J 

and a related conclusion is that a function / belongs to A4 if and only if at admits 
a representation 



n—i I n _ j— jg \ j 

/(z) = 6 a (z) • V _ J" 1 ■ 5 „_,(z) for some g k £ A k . (7.46) 

The second approach: Since the function z i-> (1 — za)^ belong to for 

n-l / _ U2 \ J 

all j, it follows that the factor 7 = — ■ g„_,(z) in (17.461) belongs to „4 n . 

3=0 \ / 

Therefore, the subspace A4 can be represented as M. — b a ■ A n ; this observation 
illustrates the second approach above. In this case, A4 with lifted norm (|7.21l) is 
contractively included in A n . 

The third approach: We first find and satisfying (|6.22p for A and C 
as in (|7.42l) and hence, with & n ,k,C,A an d <& n ,k+i,C,A defined as in (|7.44[) . We are 
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seeking Bk and Dk subject to 

l Bfe l = « 5T = 7T^T 4r ' ^ 

©n,fe+l,C,4 *&n,k,C,A ^n,k , &n.k,C,A'<->n.k+l,C,A 

Wf.^.^m. =„„,,. W!|Wd, (7 . 48) 

^ = - 5 ^-H 2 )» . (7.49) 

Note that the second equalities in (17.481) and (|7.49p follow from (|7.44|) ; alternatively 
one can use (|7.45[) combined with (j2.7[) . The essentially unique choice of Bk and 
D k satisfying (|7~S7|) - ([7~igj) is the following: 



Bk = « . a-H 2 )* = H . Maw , (7 50) 

Now we arrive at the formula for G nj fe: 

©«,*(*) = ^~fc£>* + zCi?„, fc+ i(zA)B fe (7.51) 



®n,k+l za(l-\a\ 2 ) n R n , k+1 (za) 



' (\a\ 2 <B n ,k+i - (1 - H 2 )™2:ai?„,fc + i(za)) 



fe+1 



2\n 



I I y/lln,k&n,k®n,k+l 

(1-I«| 2 )" 
|«|\/AtTi,fe'J5n,fc®n,fe+l 



(|a| 2 i?„,/c+i(|a| 2 ) ~ zaRn,k+\{zoL)) 
(RnAN 2 ) - Rn,k{za)) , (7.52) 



where we used (|7.45|) for the fourth equality and (|2.7[) for the fifth. For k = 0, the 
latter formula gives 

H - l<rR n 

e -° (z) = 1 1 y h, • ( (i - |a|2r " 

1 f l _fl-\a\ 



Oi\y/® n ,l V V 1 - Za 



n-1 /, i | 2 \ 3 



Vl-(1-M 2 ) ,-„ 

For k > 1, we use the identity 



^) = ,yp-K ) (7 53) 

a |2)n f-J I 1-zSy 



4(5 



J. A. BALL AND V. BOLOTNIKOV 



and formula (|2.5j) to get 

n r . 

i?„, fc (M 2 )-i?„, fc (za) = ]T( £ +V N 



n 

= -a& Q ( z )E ( 
£=i 


4 

£+fc-2 
£-1 


= 1 

n+i-e 

) E 

r=l 


n 

= -a6 a (z) ( 
1=1 


£+fc-2 
£-1 


n-j? 
j=0 y 


n-l 

= -ab a (z) ^ 


/n-j 


i+k-2 \ 

e-i ) 


a&Q,(z) 


ri — 1 / 
3=0 \ 


n—j—l 

E ( 


(1 - |a| 2 )" 



(1 - |a| 2 )»+ 1 -^ (1 - za) n+1 - 1 
1 



|2Wl - ,7v^+i-^-f 



1 



(1 - |a| 2 )"+ 1 -i-^ / (1 - 



M 2 ) 



1 — za 



which being substituted into (|7.52[) gives 



= -, i / #=i ; E f E m a - h 2 /) • C-r^" 

The fourth approach: Construction of the Bergman inner function associated 
with the wandering subspace for M. amounts to the construction of 0o already 
derived above in (|7.53j) . It is easily checked that the formula (|7.53[) agrees with [15] 
formula (6) page 125] (with n = 2 in (|7.53jl ). and with the formula for G a (z) in [19} 
page 58] (with a = n — 2), up to a constant multiplicative factor. The derivation 
in both [15] and [19] is based on the work of Hedenmalm [20] whereby Bergman 
inner functions are produced as solutions of an appropriate extremal problem. We 
plan to discuss how our state-space methods can be used to solve such extremal 
problems directly in a future publication. 



8. Connections with time-varying linear systems theory 

Here we make more explicit the connections of a Bergman-inner family {Q n ,k}k>o 
and the associated realization formulas (see Step 3 in the algorithm in Theorem l7.8[) 
with the theory of conservative/dissipative time- varying linear systems as presented 
e.g. in [3] and [15] . 

We suppose that we are given an exactly n-observable pair (C, A) (A s C(X), 
C 6 £(X,y)). Then Proposition 15.71 assures us that all the gramians <8 n ,k,C,A 
(k = 0, 1, 2, ... ) are strictly positive definite. Let us introduce Hilbert spaces Xk 
(k = 0, 1, 2, . . . ) with Xk — X for all k but with Xk given the inner product induced 
by the fc-th shifted gramian: 

(x,x') Xk := (<& n ,k,c,A x > x ')x- 
Similarly we set 3^fc = y with 34-inner product given by 
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In addition suppose that operators Bk G C(JAk,X) and Dk & C{U k ,y) are con- 
structed as in Step 2 of the algorithm in Theorem 17.81 We then let Ufc be the 
colligation matrix 



V k = 



C fc 

k+n r 
k+1 1 



B fc 


I 



k+n A 
k+T A 

c 



B k 



(k+n\ n 

rr 1 )^ 



(k+n-l\ t 





Xk 




%k+i 






— > 


_ y _ 



(8.1) 



(8.2) 



Here A& = fry^4 with A considered as acting from X k into Xk+i an d similarly 
Ck = C but considered as acting from Xk into y while Bk is considered as acting 
from Uk into Xk- The result of the construction in Theorem l7.8l is that each Uk '■— 
[ C dI ] * s umtar y horn Xk © Uk to Xk+i © 34- One can then study the associated 
conservative discrete-time time- varying linear system 

x(k+l) = A k x(k) + B k u(k) 
y{k) = C k x(k)+D k u(k). 

If we specify an initial condition x(0) = 0, then the input string {u(fc)}fc>o re- 
cursively uniquely determined via the system equations (|8.2|) an output string 
{y{k)}k>o, thereby defining a linear operator T{ Ufc } : {u(k)} k >o H> {y(k)} k >o, the 
input-output operator for the linear system £{u fc } fc>0 . Such systems (with k run- 
ning over all of Z rather than just Z+) are studied in [3] from the point of view 
of model theory for a sequence of contraction operators Tk = A* k : X k +x — > Xk and 
associated time- varying unitary dilation theory and Lax-Phillips scattering theory. 
However the time- varying system associated with the realization of our Bergman- 
inner family of transfer functions is not given by {Uk} but rather by {\Jk}k>0y 
namely: 

x{k + l) = A k x(k) + B fe u(fc) 
y(k) = Ckx(k)+-D k u(k). 



(8.3) 



where XJk = [ ] i s ^ ne tweak of Uk given by (|8.1|) . The input-output operator 

T{u k y k>0 associated with the system (|8.3|) gives rise to an input-output operator 
T {u,} fc > : £ 2 ({Uk}k>o) -> £ 2 (y) with image space T {Uk}k ^ (£ 2 {{U k } k >o)) exactly 
equal the set of all Taylor-coefficient strings {y(k)}k>o for which the associated 
/^-transform y(z) = J2kLo v(k)z k is in the prescribed closed S^-invariant subspace 
M of A n (y) (assuming that M is chosen so that it is generated by its wandering 
subspace). Moreover ?{u} fc>0 is isometric if the output string {y(fc)}fc>o is given 
the Bergman-space norm of its Z-transform | X^^Lo y(k) z \\A„(y)- It is amusing 
that this tweak U k i— > Ufc as in (|8.1|) of a standard conservative time- varying linear 
system ^{u k } k>0 is the precise object required to get a realization theory for a 
Bergman-inner family {0fc}fc>o- 
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